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ABSTRACT 

External  burning  increases  base  pressure.   The  character- 
istic and  linearized  equations  with  heat  addition  are  discuss 
ed,  and  the  condition  for  maintaining  irrotational ity  with 
heat  addition  is  derived.   A  transition  model,  which  extends 
the  Crocco-Lees  base  flow  model,  and  a  spin  model  are  dis- 
cussed.  Computed  results  are  compared  with  experiment.   The 
optimum  location  for  supersonic  external  burning  is  given. 


TABLE  OF  CONTENTS 


I.  INTRODUCTION 12 

II.  THEORY 18 

III.  EXPERIMENTAL  INVESTIGATION 3  7 

IV.  NUMERICAL  INVESTIGATION 48 

V.  DISCUSSION  OF  RESULTS 54 

VI  .   CONCLUSIONS 73 

APPENDIX  A:   DERIVATION  OF  THE  CHARACTERISTIC  EQUATIONS 

WITH  HEAT  ADDITION 7  5 

APPENDIX  B:   DERIVATION  OF  CONDITION  FOR  IRROTATION- 

ALITY  WITH  HEAT  ADDITION 79 

APPENDIX  C:   BRIEF  DEVELOPMENT  OF  CROCCO-LEES  BASE 

FLOW  FORMULATION 8  3 

APPENDIX  D:   CALCULATION  OF  BOUNDARY  LAYER  MASS  FLUX 

AT  THE  BASE 8  7 

APPENDIX  E:   MASS  FLUX  CALCULATION  ALONG  THE  SHEAR 

LAYER 90 

APPENDIX  F:   CALCULATION  OF  BOUNDARY  LAYER  THICKNESS 

AND  MIXING  COEFFICIENT 

APPENDIX  G:   DERIVATION  OF  LINEARIZED  STEADY  FLOW 

EQUATION  WITH  HEAT  ADDITION 96 

APPENDIX  H:   DETERMINATION  OF  THE  MACH  ANGLE 

APPENDIX  I:   SPINNING  CYLINDER  DESIGN 101 

APPENDIX  J:   COMPUTER  PROGRAM 10  3 

COMPUTER  PROGRAM 10  7 

LIST  OF  REFERENCES 1  ; .; 

INITIAL  DISTRIBUTION  LIST 128 


LIST  OF  TABLES 


I.  Base  Pressure  Ratio  and  Specific  Impulse  with 
Simulated  Heat  Addition 46 

II.  Comparison  of  Calculated  Base  Pressure  Ratio  with 
External  Heat  Addition,  Two  Dimensional  Planar 

Flow  and  Mach  Number  of  2.3 60 


LIST  OF  FIGURES 

1.  External  Burning  Projectile 13 

2.  Base  Pressure  Ratio  and  Flow  Angle 15 

3.  Natural  and  Characteristic  Coordinate  Systems 19 

4.  Conceptual  Diagram  Depicting  Shear  Layer 
Pressures,  P., ,  Which  Straddle  the  Critical 

Point 7. 2  6 

5 .  Transition  Model 2  7 

6.  Transition  of  Shear  Layer  Mixing  Coefficient 31 

7.  Momentum  Control  Volume  Depicting  Shear  Layer 

Mixing 32 

8.  Sensitivity  of  Base  Pressure  Ratio 34 

9.  Computed  Values  of  Pu/p^  vs.  Reynolds  Number, 

Mach  2.0,  No  Meat  Addition 35 

10.  Heat  Zone  Simulation  with  Nozzle  Wall 38 

11.  Photograph  of  Free  Jet  Experimental  Facility 39 

12.  Schematic  Diagram  of  Free  Jet  Experimental 

Facility 40 

13.  Coaxial  Epoxy  Nozzle.   Insert  Showing  Spacers 42 

14.  Schlieren  Photographs  at  Mach  2.0  Showing 

Uniform  Flow  and  Flow  with  Compression 44 

15.  Characteristic  Mesh 49 

16.  Integration  in  the  F  -  w   Plane 51 

"  c 

17.  Computed  Values  of  Pi/p^j  Transition  Points  and 
Critical  Point 55 

18.  Experimental  Values  of  p,  /pOT  vs.  Reynolds 

Number  Without  Heat  Addition,  from  kef.  32 56 

19.  Experimental  Values  of  p,  /p^  vs.  Reynolds  Number 
Without  Heat  Addition.  .  .  .3 .  .  ? 


20.  Comparison  of  Experimental  and  Computed  Values 
of  Base  Pressure  Ratio  vs.  Reynolds  Number  at 
Mach  2.0  Without  and  With  Simulated  Heat 

Addition 62 

21.  Comparison  of  Experimental  [Ref.  26]  and  Computed 
Values  of  Base  Pressure  Ratio  vs.  Reynolds 

Number  With  Spin 64 

22.  Computed  Values  of  Base  Pressure  Ratio  at  a  Mach 
Number  of  2.0  With  External  Burning 66 

23.  Comparison  of  Calculated  Base  Flow  With  And 
Without  Heat  Addition  at  Mach  2.0 68 

24.  Schematic  Diagram  of  Resulting  Base  Flow  When  an 
External  Shock  Wave  Impinges  on  the  Shear  Layer...  71 

II.   Disassembled  View  of  Spinning  Cylinder 102 

1 J .   Simplified  Flow  Diagram  for  External  Burning 

Assisted  Projectile  (EBAP)  Computer  Program 105 


NOTATION 

English 

a    speed  of  sound,  feet/second 

Cr    skin  friction  coefficient,  dimensi onless 

c    specific  heat  at  constant  pressure, 
"    foot-pounds/slug- degree  Rankine 

F    dimensionless  ratio  of  various  boundary 
layer  thicknesses 

h  specific  enthalpy,  foot-pounds/slug 

k  shear  layer  mixing  coefficient,  dimensionless 

L  body  length,  feet 

M  Mach  number,  dimensionless 

m  mass  flux,  slugs/second-square  foot 

n    normal  coordinate,  feet;  also 
velocity  profile  integer 

p     static  pressure,  pounds/square- foot 

Q    heat  added  per  unit  mass,  foot-pounds/slug; 

also  dimensionless  heat  addition  rate  of  Eqn.  (12) 

q    heat  addition  rate,  foot -pounds/slug -second 

R    universal  gas  constant,  foot -pounds/slug-degrec 
Rankine 

Re  Reynolds  number,  dimensionless 

r  radial  coordinate,  feet 

S  specific  entropy,  foot-pounds/slug-degree  Rankine 

s  streamline  coordinate,  feet 

T  static  temperature,  degrees  Rankine 

U  velocity  vector,  feet /second 

U,u  speed,  feet/second 


w    speed  along  streamline,  feet/second;  also 
nondimensionalized  by  stagnation  speed 

of  sound,  a 

'   s 

w/r   dimensionless  heat  zone  length 

z/r   dimensionless  position  on  body  or  shear  layer, 
measured  from  base 

Greek 

a  flow  field  type,  dimensionless 

3  (M2   -  l)h 

Y  ratio  of  specific  heats,  dimensionless 

6  boundary  layer  thickness,  feet 

6*  boundary  layer  displacement  thickness,  feet 

6**  boundary  layer  momentum  thickness,  feet 

c  dimensionless  small  quantity,  0  <  e  <<  1 

n     left  running  characteristic  coordinate,  feet; 
also  dimensionless  heating  parameter  of  Ref.  7 

6     flow  angle,  radians 

k     dimensionless  ratio  of  average  internal  to 
external  speed 

u  Mach  angle,  radians 

v  Prandtl -Meyer  angle,  radians 

£  right  running  characteristic  coordinate,  feet 

p  density,  slugs/cubic  foot 

cf>  perturbation  velocity  potential  ,  square  feet/second 

Q  curl  U,  radians/second 

a)  angular  velocity,  radians/second 

Superscripts 

(  )'  dimensionless  perturbation  quantity 
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Subscripts 

1  average  condition 

00  freestream  condition 

b  base  condition 

e  condition  at  edge  of  shear  layer 

L  laminar 

o,s  total  or  stagnation  condition 

T  turbulent 

tr  laminar  transition  point 

trt  turbulent  transition  point 
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I.   INTRODUCTION 

Projectiles  in  supersonic  flight  experience  drag.   One 
source  of  drag  is  due  to  the  pressure  acting  on  the  base 
(rear  surface)  of  the  projectile,  which,  under  normal  cir- 
cumstances, is  below  ambient  pressure.   The  base  drag  force 
is  equal  to  the  base  surface  area  normal  to  the  direction 
of  flight  multiplied  by  the  mean  difference  between  the 
ambient  and  base  pressures. 

External  burning  has  been  studied  by  many  authors  [Ref. 
1]  to  determine  its  ability  to  reduce  the  drag  of  a  body  or 
even  to  produce  thrust.   The  basic  objective  of  external 
burning  in  a  supersonic  flow,  Fig.  1,  is  to  create  compres- 
sion waves  which  impinge  on  the  base  flow  region  and  increase 
the  static  pressure  within  the  recirculation  zone,  thereby 
increasing  the  pressure  acting  on  the  base  of  the  body  to 
reduce  base  drag  or  produce  thrust. 

Other  methods  of  reducing  base  drag  or  producing  thrust 

are  available.   The  rocket  assisted  projectile  provides  thrust 

from  combustion  of  self -contained  fuel  and  oxidizer.   With 

external  burning,  additional  fuel  can  replace  some  of  this 

oxidizer  since  oxygen  is  used  from  the  atmosphere  for  burnii 

the  excess  fuel.   The  ram  jet  assisted  projectile^ delivers 

excellent  performance  but  docs  require  air  inlets  and  has 

additional  hardware  costs.   The  method  of  base  burning  offers 

advantages  similar  to  external  burning,  but  has  not  proven 

its  ability  to  create  thrust. 
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The  external  burning  assisted  projectile  provides  a  509o 
increase  in  range  when  fuel  rich  propellant  is  substituted 
for  one  third  of  the  conventional  projectile  payload.   In 
addition  to  an  increase  in  range,  external  burning  can  be 
used  to  decrease  the  projectile  flight  time  to  the  target, 
increase  its  kinetic  energy  at  the  target  and  provide  a 
flatter  trajectory. 

The  base  pressure  ratio,  p, /p^,  is  a  measure  of  the  ef- 
fectiveness of  the  external  burning  process.   As  shown  in 
Fig.  2,  the  behavior  and  effect  of  this  ratio  can  be  related 
to  the  flow  angle,  6,  of  the  boundary  layer  streamline  at 
the  base. 

For  purposes  of  analysis,  the  external  burning  process 
may  be  described  thermodynamically  as  a  heat  addition  term 
in  the  energy  equation,  bearing  in  mind  that  numerous  other 
factors  must  be  considered  with  actual  combustion,  such  as 
the  effect  of  pressure  on  combustion  rate.   The  linearized 
problem  involving  two  dimensional  planar  uniform  flow  with 
unit  heat  addition  has  been  solved  by  Tsienand  Beilock  in 
closed  form  [Ref.  2,3],   In  supersonic  flow  it  was  shown 
that  a  pressure  perturbation,  due  to  the  heat  source,  is 
transmitted  along  a  characteristic.   Further  study  of  the 
linearized  flow  equations  by  Rues  has  shown  that  there  is  an 
equivalence  between  heat,  force  and  mass  sources  [Ref.  4]. 
These  unit  sources  were  given  an  integral  formulation  by 
Fuhs  to  determine  the  flow  field  perturbation  resulting  from 
a  source  zone  [Ref.  5],   It.  is  the  study  of  these  Lineari 
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results  which  suggested  that  for  experimental  investigation 
the  external  burning  process  may  be  simulated  by  supersonic 
nozzle  wall  contours  which  generate  identical  characteristics 

Heat  addition  in  a  two  dimensional  planar  supersonic  flow 
was  investigated  by  Strahle,  using  one  dimensional  and  lin- 
earized heat  addition  theory,  and  shown  to  reduce  base  drag 
or  produce  thrust  [Ref.  6,7].   In  this  study,  the  base  flow 
was  described  by  the  Crocco-Lecs  [Ref.  8]  model,  and  the 
resulting  base  pressure  ratio  was  determined  for  heat  addi- 
tion zones  of  different  intensity,  location  and  extent.   Ex- 
perimental investigation  [Ref.  9-11]  has  demonstrated  the 
feasibility  of  external  burning  and  shown  that  it  can  in- 
crease base  pressure.   Heat  addition  in  an  axisymmetric 
supersonic  flow  has  been  studied  by  Hosack  and  O'Leary  using 
the  method  of  characteristics  and  the  Korst  base  flow  model 
[Ref.  11] . 

The  object  of  the  present  investigation  was  to  study 
base  pressure  in  two  dimensional  planar  and  axisymmetric 
supersonic  flow  with  external  heat  addition  using  the  method 
of  characteristics  and  the  Crocco-Lees  base  flow  model.   The 
effects  of  varying  Reynolds  number,  including  boundary  layer 
transition,  and  projectile  spin  have  been  included.   A  con- 
current experimental  investigation  studied  the  effect  on 
base  pressure  of  various  axisymmetric  nozzle  contours  which 
simulated  the  external  burning  process. 

It  should  be  mentioned  that  several  combinations  o\~    ex- 
ternal heat  addition  analysis  and  base  flow  models  have  not 
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been  investigated  theoretically.   These  include  two  dimen- 
sional planar  method  of  characteristics  and  linearized  theory 
with  the  Chapman-Kors t  base  flow  model,  and  linearized  axi - 
symmetric  flow  with  both  the  Chapman-Korst  and  Crocco-Lees 
base  flow  models.   The  linearized  axisymmetric  flow  problem 
will  be  discussed  later. 
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II.   THEORY 

The  theory  employed  for  this  investigation  must  describe 
the  steady  flow  in  two  entirely  different  flow  regimes, 
Fig.  1.   The  external  flow  is  supersonic  and  can  be  accurately 
treated  by  the  method  of  characteristics.   The  internal  vis- 
cous flow,  or  recirculation  zone,  is  dissipative  and  diffi- 
cult to  describe.   As  in  the  two  dimensional  case,  it  was 
assumed  that  the  Crocco-Lees  model  [Ref.  1,12-15]  adequately 
describes  the  axisymmetric  internal  flow.   Boundary  layer 
type  (laminar  or  turbulent) ,  thickness  and  transition  were 
included  as  a  function  of  Reynolds  numer  since  these  effects 
profoundly  influence  the  base  flow  behavior. 

The  characteristic  equations  have  been  derived  for  flow 
with  heat  addition  [Ref.  16,17].   In  vector  notation  the 
steady  equations  for  flow  without  viscosity,  heat  conduction 
or  body  forces  are 

V-(pU)  =  0  (1) 

pU-VU  +  Vp  =  0  (2) 

U-V(h  +  yJ-U)  -  q                                  (3) 

p  =  PRT  (4) 

where  p  is  density,  U  is  velocity,  p  is  static  pressure,  h  is 
enthalpy,  q  is  the  time  rate  of  heat  addition  per  unit  mass, 
R  is  the  universal  gas  constant  and  T  is  static  temperature. 
Using  a  natural  coordinate  system  [Ref.  18],  Fig.  3,  the 
characteristic  equations  were  derived  from  the  above  equations 

18 
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in  the  following  convenient  form  (Appendix  A) 

a(v-O)    ,  asinO  rbq  ,   .      cosyr..,8S  8h(K  rc, 

K    s  -  =  (  7 -  „„  rp  )  siny  -  -ttt— (Itt-  -  — ^— J  (5) 

3n     v   r/r,  Uc  i  *                      \il    K    dn  8ny  v  J 

b  p 

8(v+6)    ,  asin0      bq  «,         cosy,,r3S      0.    r,> 

8£      v   r/r,    Uc  T  ^  U^  v  3n     8n^    v  J 

*  b     p 

where  9,  v  and  y   are  respectively  the  flow,  Prandtl -Meyer 
and  Mach  angles,  U  is  speed,  S  is  specific  entropy  and  h^ 
is  total  enthalpy,  hu  =  h  +  ^U2 .   The  left  and  right  running 
characteristics  are  r\    and  E,    respectively.   The  natural  co- 
ordinates, s  and  n,  are  along  and  normal  to  the  streamline. 
For  two  dimensional  planar  flow,  a  is  0,  and  a  is  1  for 
axisymmetric  flow.   The  characteristic  equations,  (5)  and 
(6) ,  have  been  multiplied  by  the  characteristic  length  of 
base  radius,  r,  ,  to  obtain  a  dimensionless  form.   Therefore, 
in  contrast  to  the  dimensional  coordinates  n ,  E,    and  n  used 
in  Fig.  3,  the  coordinates  n,  5  and  n  in  Eqns .  (5)  and  (6) 
are  understood  to  be  dimensionless.   Also,  the  specific  heat 
at  constant  pressure,  c  ,  was  assumed  constant  in  the  enthalpy 

relation  h  =  c  T. 
P 

The  characteristic  computations  with  Eqns.  (5)  and  (6) 

1   9S    9h0 
were  simplified  by  neglecting  the  term  ff(T^—  -  — ^—  )/r,  , 

which  is  the  vorticity.   One  may  assume  the  external  flow  to 

be  irrotational  before  heat  addition  and  then  determine  the 

restrictions  on  heat  addition  which  preserve  irrotationali ty . 

Crocco's  theorem  for  steady  flow  may  be  written  (Appendix  B) 

TVS  +  UXTT  =  Vh  (7) 

where  H   =  VXU.   Taking  the  curl  of  Eqn.  (7) 
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U'Vft  =  Q-VU  -  ftV-U  +  VTXVS.  (8) 

Assume  fi  =  0  before  and  during  heat  addition.   Then,  from 
Eqn .  (8),  the  condition  which  preserves  irrotati onali ty  is 

VTXVS  =0.  (9) 

If  one  redefines  the  heat  addition  and  has  uniform  upstream 
conditions  at  infinity  [Ref.  19],  then 

V(hQ  -  Q)  =  0  (10) 

where  Q  is  the  heat  per  unit  mass  which  has  been  supplied  to 
the  fluid  element.   Using  Eqn.  (7),  with  Q   -  0,  Eqns .  (9) 
and  (10)  yield  the  following  result 

VQXV(yj2)  =  0  (11) 

which  is  the  condition  preserving  irrotationality ;  that  is, 
the  gradients  of  kinetic  energy  and  heat  added  must  be  par- 
allel.  If  Q  is  a  small  perturbation,  (Q/h  )  =  c,  then  Eqn. 
(11)  is  of  order  e2  and  the  irrotationality  condition  is 
automatically  satisfied  to  order  e.   To  neglect  the  vorticity 
term  in  Eqns.  (5)  and  (6),  one  must  insure 

OT-T  =  rbl7cT  "  *  «  1.  C12) 

P        P 

In  the  computations  discussed  in  Section  V,  the  quantity 
of  Eqn.  (12),  (r,q)/(Uc  T) ,  was  assigned  a  constant  value 
independent  of  local  pressure,  within  the  heat  zone  and  was 
taken  as  zero  outside  the  heat  zone.   Thus,  within  the  heat 
zone,  VQ  was  approximately  parallel  to  the  gradient  of  kinetic 
energy  and  the  irrotationality  condition  of  Eqn.  (]I)  was 
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satisfied.   This  condition  was  not  satisfied  along  the  inner 
and  outer  streamlines  of  the  heat  zone,  since  a  jump  in  Q 
existed  normal  to  the  streamlines.   However,  by  assigning  a 
value  for  Eqn.  (12),  (r,q)/(Uc  T)  <  0.2  <<  1,  the  irrota- 
tionality  condition  was  satisfied  to  order  epsilon  along  the 
streamline  boundaries  of  the  heat  zone. 

Equations  (1)  through  (4)  may  be  combined  to  give 

VU  -  ^.VOsU-U)  =  ^_  (13) 

P 

where  a  is  the  speed  of  sound.   Assuming  small  perturbation, 
Eqn.  (13)  may  be  linearized  (Appendix  G) 

v2*  -  K  0  ■  q"  (") 

where  Mto  is  the  Mach  number  at  infinity,  q'  is  the  heat  ad- 
dition, q  =  h  q1,  <J>  is  the  perturbation  velocity  potential, 
U  =  lU   +  tJ!  ,  U'  =  V<j>  and  T  =  T  fl+T1).   In  Cartesian  co- 

OO  *  T  OO  v  J 

ordinates ,  Eqn.  (14)  may  be  written 

-    32cf>        +   <f>        +    (f>        =   q'  (15) 

xx  yy        Yzz         H  K      J 

where  $2  =  M2  -  1.   For  supersonic  flow,  the  solution  of 

OO  l  ' 

Eqn.  (15)  is  found  by  integrating  the  sources  in  the  fore- 
cone  volume,  Vr  [Ref.  20], 

.       1    I q'dxidyidzi nfi. 

*  "  "  2?  \f       ((x-Xl)2  -  B2((y-yi)2  -    (z-z,)^-  UU 
vf 

Although  the  linearized  theory  was  not  pursued  further,  it 
should  be  noted  that  this  theory  could  be  used  instead  of 
the  method  of  characteristics  to  establish  the  effect  of 
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external  burning  at  the  base  flow  shear  layer.   Studies  of 
linearized  supersonic  f]ow  through  ducts  have  provided  solu- 
tions for  an  annular  source  distribution  [Ref.  21-23],  which 
could  be  used  in  conjunction  with  centerline  sources  and 
sinks  [Ref.  18]  to  describe  the  axisymmetric  external  flow 
with  burning. 

The  two  dimensional  planar  Crocco-J.ecs  theory  [Ref.  8] 
considers  the  base  flow  as  a  separated  boundary  layer  merg- 
ing into  a  free  jet  profile  and  derives  integral  equations 
to  describe  the  mean  values  of  the  flow  quantities  (Appendix 
C) .   From  the  integrals  of  the  conservation  equations,  a 
single  ordinary  differential  equation  is  derived.   This  or- 
dinary differential  equation  has  a  saddle  point  singularity; 
existence  of  the  singularity  is  a  key  feature  of  the  Crocco- 
Lees  model.   The  behavior  of  this  equation  is  partly  deter- 
mined by  an  empirical  F  -  k  relation  and  by  the  nature  of 
the  shear  layer  mixing.   k  is  the  ratio  of  average  internal 
base  flow  velocity  to  the  external  velocity,  u./u  .   F  is  a 
ratio  formed  by  using  the  various  boundary  layer  thicknesses. 
The  turbulent  shear  layer  mixing  coefficient,  k™ ,  was  as- 
signed the  constant  value,  0.03,  by  Crocco-Lees.   The  laminar 
value  used  for  this  investigation  was  0.0033,  about  1/9  the 
turbulent  value. 

F  is  conceptually  a  profile  form  factor;  when  the  veloc- 
ity profile  changes  shape,  the  value  of  F  changes.   As  dis- 
cussed in  reference  8,  the  F  -  k  functional  relationship  is 
based  on  many  known  viscous  flows.   For  example  the  Blasius 
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solution  is  represented  by  a  point  in  the  F  -  K  plane. 
Falkner-Skan  boundary  layer  solutions,  which  involve  either 
favorable  or  adverse  pressure  gradients,  become  a  curve  in 
the  F. -  K  plane.   The  velocity  profiles  of  jets  under  cer- 
tain conditions  are  described  by  the  error  function;  these 
jets  become  a  point  in  the  F  -  k  plane.   Likewise  the  Howarth 
boundary  layer  solutions  for  flows  in  pressure  gradients  are 
a  curve  in  the  F  -  k  plane. 

The  F  -  k  curve  corresponds  to  viscous  flow  near  walls, 
i.e.  boundary  layers,  as  well  as  free  jets.   Furthermore, 
adverse,  zero  or  favorable  pressure  gradients  are  represented 
by  positions  on  the  F  -  k  curve.   Heat  addition  in  a  super- 
sonic flow  generates  compression  waves.   When  z/r  <  0,  the 
compression  waves  impinge  ahead  of  the  base;  the  boundary 
layer  on  the  body  and  the  shear  layer  experience  an  adverse 
pressure  gradient.   When  z/r  >  0,  the  compression  waves  im- 
pinge downstream  of  the  base  and  only  the  shear  layer  flows 
into  the  adverse  pressure  gradient  established  by  the  exter- 
nal burning.   Since  the  F  -  k  curve  incorporates  pressure 
gradients,  the  presence  of  an  adverse  pressure  gradient  due 
to  heat  addition  should  not  alter  its  validity.   As  stated 
in  reference  8,  the  F  -  k  curve  needs  additional  experimental 
verification;  this  applies  to  the  case  with  external  burning. 
Further  discussion  of  the  F  -  k  curve  may  be  found  in  refer- 
ences 6,  12,  13,  14  and  15. 

The  reformulation  of  the  Crocco-Lees  theory  for  axisym- 
metric  flow  [Ref.  1,12-14]  resulted  in  the  same  single  ordinary 
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differential  equation  which  had  been  derived  for  planar  flow. 
In  two  dimensional  planar  flow  this  equation  is  numerically 
integrated  by  starting  at  the  singularity,  or  so  called 
critical  point,  and  working  toward  the  base,  using  the 
Prandtl -Meyer  relation  to  obtain  flow  angle  for  various  Mach 
numbers.   With  external  heat  addition,  one  may  no  longer  use 
the  Prandtl -Meyer  relation  for  flow  angle  at  the  shear  layer, 
since  the  heat  addition,  described  by  either  linearized  or 
one  dimensional  flow,  introduces  a  correction  term  [Ref.  6,7] 

This  procedure  cannot  be  used  for  the  axisymmetric  prob- 
lem without  heat  addition  due  to  the  1/r  term  arising  in  the 
continuity  equation;  thus,  one  must  start  at  the  base  and 
integrate  toward  the  critical  point.   As  was  anticipated, 
attempts  to  numerically  integrate  through  this  saddle  point 
singularity  were  unsuccessful;  however,  by  repeated  integra- 
tion one  may  determine  two  Mach  numbers  on  the  downstream 
side  of  the  expansion  fam  which  straddle  the  saddle  point 
singularity,  thus  establishing  the  base  pressure  within  the 
limits  given  by  the  two  Mach  numbers,  Fig.  4.   This  technique 
was  used  with  Eqns .  (5)  and  (6)  for  both  axisymmetric  and 
two  dimensional  planar  flow. 

To  incorporate  the  effects  of  various  boundary  layers  on 

the  body  (laminar,  transitional,  turbulent)  and  shear  layer 

transition,  it  was  necessary  to  formulate  a  transition  model, 

Fig.  5.   While  this  model  treats  the  boundary  layer  mass  flux 

differently,  it  is  not  a  departure  from  the  existing  Crocco- 

Lees  theory.   In  fact,  this  transition  model  extends  the 

Crocco-Lces  base  flow  model. 
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•Since  the  transition  from  laminar  to  turbulent  flow  is 
more  or  less  probabilistic  in  nature  and,  depending  on  the 
Reynolds  number,  may  occur  on  the  body  or  along  the  shear 
layer,  it  was  assumed  that  the  transition  process  could  be 
described  by  an  integrated  normal  distribution  and  two 
transition  Reynolds  numbers.   The  transition  Reynolds  num- 
bers .identified  the  start  and  end  of  the  transition  process, 
and  the  integrated  normal  distribution  ranged  from  0.02  to 
0.98  in  the  span  between  the  two  specified  transition 
Reynolds  numbers.   This  assumption  replaces  the  previously 
assumed  transition  curves  [Ref.  8]. 

The  transition  model  must  provide  several  parameters 
which  are  required  in  the  Crocco-Lees  base  flow  formulation. 
It  must  establish  the  values  for  the  mixing  coefficient,  k, 
along  the  shear  layer  and  the  boundary  layer  mass  flux  at 
the  base  to  permit  determination  of  F. 

The  boundary  layer  was  assumed  to  be  at  constant  static 
pressure  and  total  enthalpy,  and  the  velocity  profile  was 
assumed  to  have  a  1/n  dependence,  where  n  is  an  integer, 
1  <  n  <  7.    This  formulation  may  be  integrated  (Appendix  D) 
for  both  planar  and  axisymmetric  flows  to  give  the  boundary 
layer  mass  flux  at  the  base  as  an  explicit  function  of  veloc- 
ity profile,  1/n,  base  radius,  r,  ,  and  boundary  layer  thick- 
ness, 6.   If  the  flow  over  the  body  is  completely  laminar, 
6  is  given  by  the  compressible  Blasius  solution  [Ref.  8], 
6L  =  (5.0  +  1. 2(Y-l)M2)h/(Rch)'5,  and  the  value  of  n  =  1  is  used 
for  the  velocity  profile  [Ref.  24].   If  the  boundary  layer  is 
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completely  turbulent,  the  n  is  taken  as  7  and  6  is  computed 

from  the  compressible  one-seventh  power  solution  [Ref  8], 

&      =    (0.037L/(keb)1/5)  (Cf  /Cf  )  (6/6**)  .   The  term  (CfM/Cf.) 

Mi  1 

is  a  ratio  of  skin  friction  coefficients  for  compressible 

and  incompressible  turbulent  flow  [Ref.  24].   It  was  recog- 
nized that  the  boundary  layer  thickness  computed  for  planar 
flow  will  be  different  from  the  axisymmetric  thickness. 
This  difference  will  be  small  if  6/r,  and  radial  velocity, 
v/U,  are  of  order  e.   The  dimensional  term  pv/r  from  the 
axisymmetric  continuity  equation  will  be  of  order  e2  upon 
nondimensionalizing  with  6  and  U.   Thus,  the  governing  equa- 
tions and  their  solutions  are  identical  to  order  e  in  both 
cases . 

The  first  transition  Reynolds  number,  Re   ,  is  where  the 
laminar  flow  starts  transition  to  turbulent.   The  second 

transition  Reynolds  number,  Re^  . .  is  where  the  flow  must  be 

'  '    trt ' 

completely  turbulent  (Re     =  4Re   )  .   As  one  changes  values 

of  base  Reynolds  number  from  Re„_   to  Re^  .   .  the  values  for  n 
'  tr      trt' 

and  6,  which  determine  the  boundary  layer  mass  flux  at  the 

base,  must  vary  from  1  to  7  and  6y  to  6T,  respectively.   This 

transition  was  computed  by  the  error  function,  erf,  since  it 

is  readily  available  in  FORTRAN  and  can  be  directly  related 

to  the  integral  of  the  normal  distribution.   A  numerical 

span,  S,  of  4,  which  takes  the  value  of  erf  from  -0.96  to 

0.96,  was  used  to  normalize  the  difference  (Rc^    -  Re4  ). 

v      trt  t  r J 

Lot   Retra   =   «totr   ♦   Retrt).      Then 

ERFB    =    ,,er£(_S    (R^^    .    Reb)/(Rctrt    -    Rctr)).  (17) 

29 


So',    when   Re,     >    Re 


r 


and 


6BASE    "    6L^    +    ERV    +    6T^    "    ERIV  ^18) 


nBASE    =    1(-H   +    ERV    +    Hh    ~    ERFB}  (]9) 


were  the  values  used  to  compute  the  boundary  layer  mass  flux 
at  the  base  (Appendix  F) . 

When  transition  from  laminar  to  turbulent  flow  occurred 
on  the  shear  layer,  the  value  for  the  mixing  coefficient, 
k ,  at  a  particular  point  was  computed  in  a  similar  fashion 
using  the  error  function  to  transition  from  laminar  mixing, 
k.  ,  to  turbulent  mixing,  kT ,  Fig.  6. 

The  effect  of  shear  layer  mixing  on  base  pressure  is 
shown  schematically  in  the  momentum  control  volume  of  Fig.  7 
The  axial  pressure  on  surfaces  cd  and  de ,  and  momentum  flux 
leaving  the  volume,  de ,  must  be  balanced  by  the  pressure  on 
surfaces  ab  and  be  and  momentum  flux  entering  the  volume, 
abed.   Thus,  to  maintain  this  balance,  an  increase  of  momen- 
tum flux  along  cd  due  to  mixing  requires  a  corresponding  de- 
crease of  base  pressure,  p,  ,  along  ab . 

One  additional  assumption  was  necessary  to  obtain  good 
agreement  with  experimental  data.   At  low  Reynolds  numbers 
and  completely  laminar  flow,  it  was  necessary  to  increase 
the  value  of  the  laminar  mixing  coefficient",  k.  ,  in  direct 
proportion  to  the  increase  of  the  thickness  of  the  boundary 
layer  at  the  base.   This  assumption  may  be  expressed  as 


i 


kL  =  k(Retr/Rebr  (20) 
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where  Re    is  the  transition  Reynolds  number,  Re,  is  the 

Reynolds  number  of  the  body,  based  on  ]ength  or  chord,  L, 

and  k  is  the  laminar  mixing  coefficient  when  Rc„   =  Re,  . 

b  tr     b 

Some  justification  may  be  provided  for  Eqn.  (20)  by  consider- 
ing the  Blasius  boundary  layer  solution  where  the  mass  flux, 
m,  is  given  by  3 . 28pux/ (Re) 2  and  the  Crocco-Lees  theory 

where  kpu  is  taken  equal  to  dm/dx.   If  one  forms  the  deri- 

U         —  u 

vative  of  3 . 28pux/ (Re)  2 ,  dm/dx  is  found  to  vary  as  l/(Re)2, 

thereby  suggesting  the  form  assumed  in  Eqn.'  (20).   Although 

kT  was  held  constant  at  0.03  for  all  calculations,  a  similar 

analysis  using  the  incompressible  turbulent  solution  suggests 

1/5 
that  kT  might  vary  as  l/(Re)  '  . 

It  has  been  suggested  that  the  effect  of  projectile  spin 

could  be  included  in  the  mixing  coefficient  [Ref.  15],  k, 

by  using  k     =  k(l  +  rco/U),  where  rio/U  is  the  spin  number. 
'      6   spin     y  '     j  >  i 

However,  since  k   .   always  increased  with  an  increase  in 
'        spin     J 

spin  number,  and,  since  it  has  been  determined  through  numer- 
ical experiments  with  the  Crocco-Lees  theory,  Fig.  8,  that 
increasing  k  decreases  the  base  pressure,  this  simple  hypoth- 
esis contradicts  experiment  [Ref.  25,26]  where,  in  some 
cases,  increasing  the  spin  number  increased  the  base  pressure 
In  Fig.  9  it  is  clear  that  at  certain  values  of  Reynolds  num- 
ber, Re,  ,  an  increase  of  Reynolds  number  resulted  in  an  in- 
crease in  base  pressure.   This  fact  suggests  that  the  spin 

number  could  be  used  to  shift  the  p,  /p   versus  Re,  curve, 

rb  l  °°  b 

Fig.  9,  instead  of  increasing  the  mixing  coefficient.   Sev- 
eral methods  were  attempted  at  this  point,  such  as  increasing 
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Reb  by  [Ref.  27]  (l+C1^-)2)^2  or  (l+C1^-)2),  but  the  method 
which  provided  the  best  correlation  was  to  simply  decrease 
the  transition  Reynolds  numbers.   The  spin  transition 
Reynolds  numbers  used  to  compute  p,  /p^  versus  Re,  were  form- 
ed by  dividing  the  corresponding  no  spin  transition  Reynolds 
numbers  by  (l+(— ry) )  •   Interpretation  of  base  pressure  data 
for  a  spinning  projectile  is  given  in  Section  V. 
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III.   EXPERIMENTAL  INVESTIGATION 

To  carry  out  experiments  related  to  the  five  inch  pro- 
jectile, the  mid  point  of  a  maximum  range  trajectory  was 
selected.   This  flight  condition  established  an  altitude 
of  23,000  feet  and  a  Mach  number  of  2.0  for  determining 
similarity  parameters.   For  the  five  inch  projectile,  this 
gave  a  body  Reynolds  number  of  15.7x10  ,  a  radius  Reynolds 
number  of  1.51x10  ,  and  a  turbulent  boundary  layer  thick- 
ness to  radius,  6/r,  of  0.1313.   The  spin  at  23,000  feet 
was  assumed  to  be  the  same  as  the  muzzle  spin,  250  revolu- 
tions per  second,  which  gave  a  spin  number,  rw/U,  of  0.1595. 
Using  these  nomial  values,  it  was  found  that  all  similarity 
parameters  could  be  matched  with  a  free  jet  Mach  2.0  coaxial 
nozzle  exhausting  at  room  pressure  if  the  center  cylinder 
had  a  radius  of  0.641  inches,  6  of  0.084  inches,  and  a  spin 
of  800  revolutions  per  second.   As  stated  in  the  introduc- 
tion, the  actual  combustion  process  was  simulated  with 
various  nozzle  contours  which  duplicated  the  compression 
waves  generated  by  the  external  burning. 

The  simulation  process  is  depicted  in  Fig.  10  where  the 
streamline  deflected  by  the  presence  of  a  heat  zone  can  be 
replaced  by  a  nozzle  wall  to  generate  an  identical  flow 
field. 

The  free  jet  experimental  facility,  Figs.  11  and  12,  was 
conventional  in  most  respects  and  has  been  previously  re- 
ported [Ref.  28,29].   The  use  of  cast  aluminum  filled  epoxy 
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(b)  Identical  Flow  Field  With  Compression  Generated  by 
Nozzle  Wall  to  Simulate  Heat  Zone 


Figure  10 ♦  Heat  Zone  Simulation  With  Nozzle  Wall 
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Figure  11.  Photograph  of  Free  Jet  Experimental  Facility 
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norzlcs,  Fig.  13,  was  the  only  new  hardware  innovation.   The 
technique  of  casting  employed  a  wax  plug  contoured  to  the 
desired  shape  with  an  aluminum  template.   This  plug  was 
turned  slightly  undersize  and  then  restored  to  size  using  a 
thin  coat  of  automotive  body  plastic.   After  melting  the  wax 
from  the  cast  nozzle,  the  plastic  peeled  away  from  the  epoxy 
to  leave  a  smooth  nozzle  wall.   Without  this  plastic,  the 
heat  generated  by  the  epoxy  had  softened  the  wax  and  caused 
some  irregularity  in  the  nozzle  wall. 

The  location  of  the  zone  of  combustion  to  be  simulated 
by  the  nozzle  was  estimated  using  an  empirical  formula  for 
penetration  depth  of  a  jet  normal  to  a  supersonic  free- 
stream  [Ref.  30].   For  Mach  2.0  flow,  internal  pressures 
from  300  to  500  psia  were  computed  to  provide  penetration 
to  2h   inches  with  various  orifice  diameters.   Thus  it  was 
assumed  that  external  burning  at  5  inches  from  the  projec- 
tile centerline  would  be  practical.   Assuming  a  fuel  rich 
propellant  exhaust,  ARC-168,  with  an  available  heating  value 
of  7200  BTU/lb  and  0.6  lb/sec  flow  rate,  it  was  found  using 
one  dimensional  theory  that  the  flow  should  be  deflected 
about  4h°    for  a  distance  of  15  inches,  or  three  diameters, 
downstream.   This  computation  used  the  integral  methods  of 
reference  5  which  gave  a  perturbation  velocity  normal  to  the 
flow  of  -163  feet  per  second.   The  resulting  flow  angle 
tangent,  was  -163/2058,  which  gave  a  computed  flow  angle  of 
-4.54  degrees  along  the  inner  streamline  of  the  heat  zone. 
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Figure  13. Coaxial  Epoxy  Nozzle.  Insert  Showing  Spacers 
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These  estimated  values  for  the  heat  zone  were  taken  as 
the  upper  limit  and  were  used  to  form  a  test  matrix  and  for 
subsequent  nozzle  design.   Several  different  nozzle  contours 
were  computed  by  the  method  of  characteristics  for  4^°  stream 
line  deflection  at  a  radius  equal  to  a  cylinder  diameter 
from  the  centerline  and  with  different  assumed  lengths  for 
the  combustion  zone.   The  test  matrix  included  using  spacers, 
Fig.  13  insert,  with  each  nozzle  so  that  the  first  full  com- 
pression waves  impinged  ahead  of,  at,  or  behind  the  base  of 
the  cylinder.   Pressures  were  measured  along  the  base  and 
side  of  a  stationary  centerline  cylinder,  Fig.  13,  using 
mercury  manometers.   Attempts  to  carry  out  experiments  with 
a  spinning  cylinder  were  unsuccessful  due  to  faulty  air  bear- 
ing design  (Appendix  I) .   Time  did  not  permit  correction  of 
this  difficulty. 

The  differences  in  the  flow  field  with  the  stationary 
cylinder  and  different  nozzles  can  be  seen  in  Fig.  14.   The 
photograph  of  the  uniform  flow  nozzle  in  Fig.  14a  used  a 
horizontal  knife  edge  and  the  expansion  fan  at  the  base  ap- 
pears dark.   The  base  pressure  ratio  was  0.63S.  In   Fig.  14b 
the  vertical  knife  edge  was  positioned  to  make  a  positive 
density  gradient  from  the  centerline  to  the  right  appear 
brighter.   The  flow  with  compression,  Fig.  14c,  had  an  ex- 
pansion fan  at  the  base  but  achieved  a  higher  base  pressure 
ratio  of  1.065  due  to  compression  upstream  of  the  base  as 
well  as  along  the  shear  layer.   The  compression  upstream  of 
the  base  increased  the  pressure  right  before  t lie  expansion 
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Figure  14a. 
Uniform  Flow, 
Horizontal 
Knife  Edge 


Figure  14b. 
Uniform  Flow, 
Vertical 
Knife  Edge 


Figure  14c. 
Compression  Flow, 
Horizontal 
Knife  Edge 


Figure  14.  Schlieren  Photographs  at  Mach  2.0  Showing 
Uniform  Flow  and  Flow  With  Compression 
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fan. so  that  it  was  greater  than  pro,  thus  the  flow  was  near 
ambient  pressure  after ■ passing  through  the  expansion  fan. 
Some  test  results  using  the  stationary  cylinder  are  given 
in  Table  I.   Except  for  the  nozzle  with  w/r  =  4.34,  these 
tests  indicate  that  compression  applied  upstream  of  the 
base,  z/r  =  -0.5,  is  less  efficient  at  increasing  base  pres- 
sure than  compression  applied  at  the  base,  or  slightly  down- 
stream.  The  nozzle  with  w/r  =  4.34  had  a  strong  shock  wave 
at  the  exit,  due  to  compression  waves  reflected  from  the 
cylinder  [Ref.  15],  which,  with  z/r  =  -0.5,  impinged  on  the 
base  flow  region  ahead  of  the  critical  point.   In  turn,  in- 
creased pressure  across  the  shock  wave  was  impressed  on  the 
subsonic  flow  in  the  recirculating  zone  and  gave  a  higher 
base  pressure.   The  shock  wave  observed  with  this  experimen- 
tal nozzle  would  not  be  expected  with  uniform  external 
burning . 

In  Table  I  the  experimental  base  pressure  ratios,  Pu/pTO, 
are  shown  as  the  upper  number  in  each  box  for  the  various 
test  positions,  z/r,  of  each  nozzle,  w/r.   The  cylinder 
radius  is  r,  the  axial  location  where  the  first  full  strength 
compression  wave  impinges  on  the  body  or  shear  layer  is  z 
(z=0  at  the  base)  and  w  is  the  length  of  the  heat  zone.   The 
corresponding  Reynolds  number,  Re,  ,  is  shown  below  each  value 
of  z/r.   The  lower  number  shown  in  each  box  is  the  computed 
specific  impulse  for  the  assumed  projectile  flight  condition 
at  23,000  feet.   The  uniform  flow  nozzle,  with  w/r  =  0.0, 
was  taken  as  the  reference  pressure  ratio  to  compute  the 
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0.0 
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0.60 

0.72 
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1.84 

0.96 
216 

1.055 
283 

1.074 
302 



3.91 

1.078 
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1.175 
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1.129 
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4.34 

1.20 
156 

1.14 
143 

1.085 

131 



Tabic  I.  Base  Pressure  Ratio  and  Specific  Impulse 
With  Simulated  Heat  Addition 
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base  pressure  force  created  by  the  other  nozzles.   The  fuel 
flow  rate  was  taken  proportional  to  (w/r)/(6.0),  since  6.0 
was  the  combustion  zone  length  assumed  for  the  fuel  flow 
of  0.6  lb/sec. 

Once  Mach  number  and  operating  pressure  and  temperature 
are  selected,  the  remaining  variable  controlling  Reynolds 
number  is  the  length  or  scale  of  the  apparatus.   The  free 
jet  facility  of  Figure  11  has  sufficient  mass  flow  to  per- 
mit duplication  of  Reynolds  number  at  23,000  feet  altitude 
provided  the  cylinder  has  sufficient  length.   Consequently 
the  nominal  Reynolds  number  was  chosen  as  15.7X10   for 
23,000  feet  to  match  capability  of  the  facility. 

When  the  design  of  the  spinning  cylinder  was  initiated 
it  became  apparent  that  a  long  cylinder  had  problems  in 
regard  to  shaft  whip  or  the  first  vibrational  mode  in 
bending.   A  shorter  cylinder  reduced  the  Reynolds  numbers 
to  the  values  shown  in  Table  I.   Conditions  at  a  higher 
altitude  were  simulated,  or  else  conditions  for  a  short 
projectile  at  23,000  feet  were  duplicated. 

Two  cylinders  could  have  been  designed.   Each  cylinder 
would  have  required  a  separate  family  of  nozzles.   A  choice 
was  made  to  use  a  short  cylinder  for  both  stationary  and 
spinning  tests. 

Variation  in  z/r  as  shown  in  the  sketch  with  Table  I 
changes  Reynolds  number  slightly.   The  variation  in  Re,  in 
Table  I  is  due  to  this  fact. 
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IV.   NUMERICAL  INVESTIGATION 

The  computer  program  for  this  numerical  investigation  was 
written  to  calculate  either  two-dimensional  planar  or  ax i sym- 
metric problems.   The  program  used  the  method  of  characteris- 
tics where  the  axisymmetric  term  in  Eqns .  (5)  and  (6), 
r,  (sinusinO) /r ,  was  deleted  for  two-dimensional  planar  flow. 
Similarly,  the  effect  of  radius  in  the  mass  flux  calculation, 
at  the  base  and  along  the  shear  layer,  was  deleted  or  re- 
tained as  required.   The  characteristic  mesh,  Fig.  15,  was 
computed  from  the  body  or  shear  layer  upstream  through  the 
heat  zone,  if  present,  to  a  reference  freestream  Mach  line 
(i.e.,  backwards  along  a  right  running  E,    characteristic). 
The  error  from  the  reference  condition  was  used  to  correct  v 
on  the  body,  since  G  was  given,  or  6  along  the  shear  layer 
since  v  was  computed  from  the  Crocco-Lees  theory.   Each  body 
or  shear  layer  point  was  iterated  until  the  error  at  the  up- 
stream  boundary  condition  was  less  than  10    radians.   At 
each  point  in  the  characteristic  mesh,  the  value  of  y  was 
determined  for  the  resulting  value  of  v  by  the  method  given 
in  Appendix  H.   In  retrospect,  this  method  of  characterist ics 
computation  could  be  made  more  efficient  if  one  could  devise 
a  scheme  which  works  with  arbitrary  points  along  the  body  or 
shear  layer  and  increments  the  characteristic  mesh  along  the 
reference  freestream  Mach  line.   Then  only  one  calculation 
would  be  required  along  each  characteristic,  rather  that  the 
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iterative  scheme  used  in  the  present  method.   This  approach 
was  considered,  but  it  was  not  clear  how  one  would  formulate 
the  calculation  for  the  corner  or  the  shear  layer  mass  in- 
flux (Appendix  E) . 

The  flow  field  was  calculated  along  the  body  to  the  base 
at  which  point  all  computed  values  along  the  last  character- 
istic were  stored.  These  values  were  used  to  reset  the  flow 
field  for  each  iteration  of  the  straddle  method. 

The  Crocco-Lees  theory  was  formulated  in  terms  of  w 

J  e 

instead  of  Mach  number,  where  w   is  the  velocity  at  the  edge 

of  the  shear  layer  divided  by  the  freestream  stagnation  speed 

of  sound,  w   =  U/a  .   To  start  the  straddle  method  a  value 
'   e     '  s 

of  w  was  chosen  for  the  downstream  side  of  the  expansion 
fan  which  was  too  large.   The  expansion  fan  was  divided  into 
ten  parts;  downstream  of  the  fan  the  computation  proceeded 
toward  the  critical  point.   Each  new  point  along  the  shear 
layer,  computed  by  the  Crocco-Lees  theory  (Appendixes  C,  E 
and  F) ,  required  several  iterations  along  the  characteristic 
to  the  upstream  Mach  line  boundary  condition,  where  the  flow 
angle  error  after  each  iteration  was  used  to  correct  the 
shear  layer  flow  angle. 

Integration  in  the  F  -  w   plane  followed  the  first  path 
shown  in  Fig.  16.   It  can  be  seen  that  this  path  falls  to 
the  right  of  the  path  passing  through  the    saddle  point 
singularity,  x,  and,  as  the  value  of  F  was  decreased,  the 
slope,  dF/dw  ,  became  negative  at  1'.   This  occurred  when 
the  denominator  used  to  form  dF/dw   became  positive  [Ref.  S], 
i.e.  both  the  numerator  and  denominator  of  dF/dw   were 
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negative  from  1  to  1* .   This  transition  from  positive  to 
negative  slope  with  positive  denominator  occurred  when  the 
initial  guess  for  w  was  too  large.   By  systematically  dc- 
creasing  the  initial  guess  for  w  ,  say  by  tenths,  one 
eventually  has  the  integration  path  from  2  to  2*,  which  is 
similar  to  the  previous  case,  except  that  the  negative 
slope  at  2'  is  associated  with  the  numerator,  used  to  form 
dF/dw  ,  becoming  positive.   Thus,  the  last  two  initial 
guesses  for  w  have  straddled  the  saddle  point  singularity. 
The  interval  between  these  two  values  may  be  divided,  say 
in  hundredths,  and  the  process  repeated  to  find  two  values 
of  w  which  more  closely  straddle  the  saddle  point  singularity 

The  heat  zone  was  specified  to  start  and  stop  at  certain 
values  of  axial  position,  z,  and  to  lie  between  an  inner  and 
outer  streamline.   The  quantity  r,  q/Uc  T  was  assigned  a  con- 
stant numerical  value,  between  0.15  and  0.175,  within  the 
zone,  which  gave  the  desired  - Ah°    deflection  to  the  inner 

streamline.   External  to  the  heat  zone  r,q/Uc  T  was  zero. 

bn   p 

Before  deflection,  the  inner  and  outer  streamlines  have 
values  of  approximately  2.0  and  3.0,  maximum  body  radius, 
respectively.   These  values  were  chosen  to  match  the  rate  of 
compression  increase  used  in  the  nozzle  design.   As  the  com- 
putation proceeded  downstream,  the  values  of  radius  and  flow 
angle  for  the  inner  and  outer  streamlines  were  computed  to 
insure  realistic  boundaries  for  the  heat  zone. 

When  a  sufficient  amount  of  properly  located  heat  was 
added  to  the  external  flow,  the  base  pressure  ratio  was 
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computed  to  be  greater  than  1.   It  is  interesting  to  note 
that  for  this  condition  the  method  of  characteristics  rou- 
tine generated  a  compression  fan  at  the  base,  rather  than 
the  usual  expansion  fan.   The  compression  fan  calculation 
mathematically  folds  the  flow  field  on  itself  [Ref.  31]. 
This  of  course  is  physically  impossible  and  is  replaced  in 
real  flow  by  an  oblique  shock.   For  values  of  p^/p^  =  1.5 
at  Mach  2.0,  the  loss  of  total  pressure  due  to  the  oblique 
shock  is  less  than  1%  ,  which  corresponds  to  an  entropy  in- 
crease, AS/R,  of  0.01.   This  resulted  in  an  error  of  less 
than  0.18°  in  the  turning  angle,  which  was  neglected.   There 
is  also  the  possibility  of  boundary  layer  separation  up- 
stream of  the  base  for  values  of  Pu/pra  greater  than  one 
[Ref.  7] .   This  effect  was  not  considered  in  the  calcula- 
tions.  The  computer  program  is  discussed  in  Appendix  J. 
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V.   DISCUSSION  OF  RESULTS 

The  computed  values  of  Pu/p^  versus  Re,  ,  shown  in  Fig. 

9,  were  based  on  the  apparent  transition  Reynolds  numbers 

from  the  uniform  flow  nozzle  experiments.   These  were 

Re^   =  1.95X106  and  Re^  .  =  7.8X106.   The  curves  display 
tr  trt  *  J 

the  same  general  features  and  approximate  values  found 

from  experiment  [Ref.  32-35].   It  should  be  noted  in  Fig.  9 

that  for  both  the  two  dimensional  planar  and  axisymmetric 

flow  the  maximum  value  of  p,/p   occurred  below  Re„.   and 

*d  r°°  tr 

the  minimum  value  of  pu/p   occurred  below  Re^  .  . 

1  b' *  °°  trt 

A  slightly  different  choice  of  transition  Reynolds 

numbers,  Re_   =  1.5X10   and  Re^  .  -  6.0X106,  resulted  in 
'    tr  trt  ' 

the  curve  of  p,  /p^  versus  Reynolds  number  shown  in  Fig.  17. 
One  may  compare  the  curves  of  p,  /poo  in  Fig.  9  and  Fig.  17 
with  the  experimental  curves  shown  in  Fig.  18.   While  the 
curves  are  by  no  means  identical,  they  do  agree  in  having 
the  same  general  shape.   Figure  19  shows  additional  experi- 
mental data  to  compare  with  Fig.  9.  The  locations  of  maximum  and 
minimum  computed  p,  /p^  are  controlled  by  the  choice  of  values 

for  Re„_   and  ReA  .  .   The  ratio  of  h    for  Re^  /Re.    used  in 
tr       trt  tr   trt 

these  computations  was  based  on  information  in  reference  8. 
This  ratio  is  an  approximation  which  varies  considerably  for 
different  bodies,  depending  on  their  shape  and  surface  con- 
dition.  The  maximum  value  of  p,  /pm   was  controlled  by  the 
choice  of  k,  ,  while  the  choice  of  k,„  controlled  the  second 
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maximum  value  of  Pu/p^  to  the  right  of  ^etrt-   The  slight  de- 
crease of  Pl/p^  to  the  right  of  this  second  maximum  is  in 
agreement  with  the  r/6  sensitivity  curve  of  Fig.  8,  since- 
increasing  Re  decreases  6  and  therefore  increases  r/6. 

Three  additional  curves  are  shown  in  Fig.  17;  these  are 
transition  and  critical  point  locations.   The  laminar  and 
turbulent  transition  points  are  found  simply  by  dividing 

base  Reynolds  number  into  Re,.   and  Re,  .  .   The  critical 
7  tr       trt 

point  location  was  taken  from  calculation  with  the  Crocco- 
Lees  equations  along  the  shear  layer  where  dF/dw   changed 
signs;  it  is  the  singularity  in  Fig.  16  which  was  discussed 
in  Section  IV.   In  going  from  small  to  large  values  of 
Reynolds  number,  the  shear  layer  between  the  critical  point 
and  the  base  gradually  changes  from  completely  laminar  to 
completely  turbulent  as  the  laminar  and  turbulent  transition 
points  move  upstream  to  the  base  and  then  onto  the  body. 
Where  this  occurs,  the  shear  layer  mixing  coefficient,  k, 
gradually  changes  from  k,  to  k„, .   Increasing  from  some  small 
value  of  Reynolds  number,  the  thickness  of  the  laminar 
boundary  layer  at  the  base  decreases  until  the  laminar 
transition  point  reaches  the  base  and  moves  onto  the  body. 
As  the  boundary  layer  at  the  base  transitions  to  turbulent 
its  thickness  increases.   The  thickness  continues  to  increase 
as  the  turbulent  transition  point  moves  forward  along  the 
body,  reaching  a  maximum  when  the  transition  point  is  near 
the  nose  of  the  body.   Beyond  this  point,  the  thickness 
decreases . 


In  Fig.  17,  the  value  of  P^/p^  starts  to  decrease  slight- 
ly before  reaching  Re   ,  while  the  laminar  transition  point 
is  still  on  the  shear  layer  and  the  base  boundary  layer  is 
laminar.   Although  the  boundary  layer  thickness  is  still 
decreasing  at  this  point  and  in  turn  should  provide  an  in- 
crease in  Pu/Poo*  the  value  of  k  along  part  of  the  shear 
layer  (between  the  laminar  transition  point  and  the  critical 
point)  starts  to  increase  from  k.  to  kT  and  consequently 
dominates  the  overall  change  in  P^/p^  with  Reynolds  number. 

The  mechanism  responsible  for  the  abrupt  minimum  value 

of  p,  /pw   in  Fig.  18  was  not  clear.   However,  this  minimum 

was  calculated  and  is  shown  in  Fig.  17  as  a  dashed  line. 

From  the  sensitivity  curve  in  Fig.  8,  it  appeared  that  k 

would  be  most  likely  quantity  responsible.   In  numerical 

experiments  this  abrupt  minimum  was  calculated  by  assuming 

additional  mixing.   The  additional  mixing  was  assumed  to 

have  a  normal  distribution  about  Re.  .  ,  i.e.  k  =  ^(1  + 

t  r  t '  T  K 

-x2/2      h 
(e      )/(2tt)  2)  .   The  argument,  x,  was  adjusted  so  that  the 

additional  mixing  fell  to  about  ten  percent  of  its  maximum 

value  in  2X10   either  side  of  Re,.  .  .   This  additional  mixing 

trt 

was  not  used  in  any  other  calculations  since  it  was  not 

understood  and  was  not  clearly  observed  in  the  experimental 

data  reported  later  in  this  section. 

Calculations  were  made  for  external  heat,  addition  in 

two  dimensional  planar  flow  to  compare  with  Strahlc's 

previously  reported  calculations  [Ref.  7].   The  values  of 

Re^   =  1.95X106  and  Re^    =  7.8X10   were  assumed  and  a 
t  r  trt 
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Reynolds  number  of  Re,  =  20.0X10   gave  the  same  value  of 
6/r  used  by  Strahle.   The  Mach  number  was  M  =  2.3.   The  heat 
addition  zone  was  started  at  the  same  point,  z/r  =  0.4,  and 
assigned  the  same  length,  w/r  =  3.0,  but  to  avoid  subsonic 
flow  with  the  method  of  characteristics  it  was  necessary  to 
increase  the  height  of  the  heat  zone  from  the  value  of  0-2r, 
used  by  Strahle  to  a  value  of  one  base  radius.   The  values 
of  Q  =  r,  q/Uc  T  were  reduced  accordingly  to  provide  the 
same  total  heat  addition  based  on  the  rectangular  heat  zone 
used  by  Strahle,  i.e.  Q  =  0.2n./tanu.   The  results  in  Table 
II  show  good  agreement  at  the  low  intensity  of  heat  addi- 
tion, but  only  fair  agreement  at  the  higher  intensities  with 
the  method  of  characteristics  program  showing  higher  values 
of  p^/p^.   Although  the  non  linearity  included  in  the  method 
of  characteristics  calculation  may  be  responsible  for  part 
of  the  higher  base  pressure  ratio,  the  main  cause  was  at- 
tributed to  expansion  of  the  heat  zone.   This  effect  was 
included  in  the  bounding  streamline  calculation  for  the  heat 
zone  in  the  method  of  characteristics  program  and  resulted 
in  a  greater  amount  of  total  heat  added  than  in  the  linear- 
ized theory  calculations. 


Linearized 
lheory1  J 

Method  of 
Char  act e  ri sties 

n 

0.0 
0.4 
0.8 
1  .2 

vb/vm 

0.420 
0  .  6  5  7 
0.900 
1.  158 

Q 

0.0 
0.16  7 
0.334 
0.5 

Pb/P« 

0.410 
0.669 
1.075 
1.45  +  0.0  2 

Table  II.  Comparison  of  Calculated  Base  Pressure  Ratio  with 
External  Heat  Addition,  Two  Dimensional  Planar  Flow  and  Mach 
Number  of  2.3. 
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Both  experimental  and  calculated  results  are  shown  in 
Fig.  20  for  a  Mach  number  of  2.0.   The  value  of  p./p^  = 
0.6314  at  Re,  =  10.4  was  the  average  of  nine  runs  with  the 
uniform  flow  nozzle.   The  probable  error  about  this  average 
was  ±  0.004,  with  a  maximum  observed  deviation  of  -  0.0104. 
In  comparing  the  experimental  values  of  Fig.  20  and  Table  I 
with  early  test  results  using  the  free  jet,  variations  as 
high  as  1%    in  PiVp^  were  found  for  similar  runs.   These 
variations  were  attributed  to  problems,  such  as  mechanical 
alignment  during  nozzle  installation  or  leaks  in  the  mano- 
meter tubing,  which  were  found  and  corrected  during  subse- 
quent testing. 

Early  test  data  were  rejected  and  were  not  reported 
herein.  As  experience  was  gained  with  operation  of  the 
facility,  the  1%    variations  were  narrowed  to  less  than  3%l . 

Cylinder  wall  temperatures  were  also  recorded  during 

these  nine  runs  using  a  thermocouple  imbedded  in  the  wall, 

1%  inches  from  the  base'.   These  temperatures  ranged  from 

5°F  to  35°F.   A  plot  of  Pv./pro   versus  temperature  showed 

no  correlation  scatter  in  p,  /pTO  due  to  temperature  was  less 

than  1.4%.   Since  the  stagnation  temperatures  of  the  air 

supply  were  not  recorded,  they  could  not  be  compared  with 

the  cylinder  wall  temperatures. 

The  values  of  Re,   =  1.95X106  and  Re4  .  =  7.8X106  were 

t  r  t  r  t 

assumed  for  the  computed  values  of  p,  /pm   shown  in  Fig.  20 
as  triangles.   The  experimental  values  of  p.  /p^  for  the 
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Figure  20.  Comparison  of  Experimental  and  Computed  Values 
of  Base  Pressure  Ratio  vs.  Reynolds  Number  at 
Mach  2.0  Without  and  With  Simulated  Heat  Addition 


uniform  flow  nozzle  without  heat  addition  showed  good  agree- 
ment with  the  computed  values.   Three  of  the  four  experimen- 
tal values  of  p,  /p^  with  w/r  =  3.91  compression  agreed  with 
the  computed  values.   The  observed  value  of  PiVp^  with 
z/r  =  -0.5  was  less  than  the  computed  value.   This  differ- 
ence was  attributed  to  an  adverse  pressure  gradient  along 
the  cylinder  wall  that  was  not  considered  in  the  calcula- 
tion.  In  view  of  this  result,  base  pressure  ratios  computed 
by  the  program  given  in  Appendix  J  should  be  considered  with 
caution  if  the  full  compression  due  to  external  burning  im- 
pinges ahead  of  the  base  on  the  projectile  body. 

Figure  21  shows  a  comparison  of  experimental  [Ref.  26] 
and  computed  values  of  base  pressure  ratio  resulting  from 
spin. 

In  the  Crocco-Lees  model  there  are  several  input  param- 
eters including  M  ,  kT  ,  k™ ,  Re^  ,  Re„.  .  and  F  -  k  relation- 
&   oo  >      l>   t  »    tr*    trt 

ship.   A  procedure  for  predicting  the  influence  of  spin  on 

p,  /p^  was  discussed  previously  at  the  end  of  Section  II. 

This  procedure  was  used  for  the  computed  values  in  Fig.  21. 

Choice  of  the  input  parameters  remains.   Preliminary 

calculations  were  conducted  with  M   =2.86,  Re,.   =  1.95X10 

oo        '    tr 

and  Re.  .  =  4.6X10  .   Values  of  p,/p   had  the  correct  trends 
trt  *d  F°° 

relative  to  experimental  values.   However,  the  calculated 
p,  /p^  for  Re,  =  2.6X10   and  zero  spin  was  about  0.54.   To 
decrease  the  calculated  Pi/po:),  the  Mach  number  was  arbitrar- 
ily increased  to  3.0.   For  Re,  =  2.6X10  ,  zero  spin  and  M^ 
=  3.0,  the  calculated  and  observed  p. /p^  values  agree.   In 
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retrbspcct,  additional  calculations  should  have  been  con- 
ducted with  increased  mixing,  i.e.  larger  k,,. ,  while  retain- 
ing correct  M^  =2.86.   In  fact,  spin  should  increase  k 
as  discussed  in  Section  II. 

A  better  procedure  for  calculating  the  influence  of  spin 
might  be  to  use  a  shift  in  transition  Reynolds  number  coup- 
led with  a  k     =  k  (1+ (roo/U) )  ,  as  discussed  by  Smithey. 
spin     k       k       /     j  j  >  j  ;   > 

Naber,  Caswell  and  Funs  [Ref.  15].  The  dashed  segment  of 
Ph/Poo  vs  •  ^eh  curvc  iR  Pig-  17  shows  the  influence  of  in- 
creased mixing  (larger  kT)  in  the  vicinity  of  minimum  p,  /pm. 

The  computed  values  shown  in  Fig.  21  required  four  hours 
time  on  the  IBM  360/6  7.   In  view  of  the  computing  time,  the 
procedure  suggested  in  the  preceding  paragraph  was  not 
accomplished. 

To  provide  some  estimate  of  the  optimum  location  for 
external  heat  addition,  a  fixed  fuel  flow  was  assumed  to  be 
burned  in  zones  of  different  length  and  location.   The  ini- 
tial radius  of  the  inner  streamline  was  2.0  and  the  outer 
streamline  initial  radius  was  3.0.   The  heat  zone  were 
assigned  integer  values  from  1  through  6  for  w/r.   The  zone 
with  w/r  =  1.0  was  assigned  a  value  of  Q  =  r,q/Uc  T  =  0.2. 
This  value  was  decreased  by  1/ (w/r)  for  the  other  five  zones 
The  heat  zones  were  positioned  along  the  axial  axis  starting 
with  compression  at  the  base  and  moving  downstream  in  in- 
teger values  of  radius.   The  results  of  two  of  these  compu- 
tations are  shown  as  the  two  lower  curves  in  Pig.  22.   The 
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gain  in  p,  /pm   for  m  =  0.15  lb/sec  is  seen  to  be  fairly  in- 
sensitive to  heat  zone  size  and  placement. 

Another  set  of  similar  calculations  were  made  based  on 
Q  =  0.8  at  w/r  =  1.0.   This  value  of  Q  corresponded  to  the 
design  fuel  flow  rate  of  0.6  pounds  per  second  at  23,000 
feet  altitude.   The  first  three  values  of  w/r,  3.0,2.0  or 
3.0  would  not  compute  because  the  external  flow  was  driven 
to  a  sonic  condition  by  the  base  pressure  rise  and  the  heat 
addition.   This  condition  is  analogous  of  thermal  choking; 
however,  in  two  dimensional  flow  both  streamtube  area  change 
and  heat  addition  drive  the  flow  toward  sonic  conditions. 
The  last  three  values  of  w/r  =  4.0,  5.0  and  6.0,  did  not 
experience  subsonic  conditions  in  the  external  flow  and 
their  values  of  base  pressure  ratio  are  shown  as  the  upper 
three  curves  in  Fig.  22.   The  high  value  of  p,  /p^  =  1.91 
for  the  w/r  =  4.0  curve  was  not  expected.   This  represents 
a  factor  of  six  increase  in  base  pressure  over  the  two 
bottom  curves  in  Fig.  22  with  a  factor  of  four  increase  in 
fuel  flow.   The  combustion  for  this  flight  condition,  Fig. 
23,  was  assumed  to  start  at  (z/r)  =  -1.2  and  the  initial 
compression  was  on  the  shear  layer,  downstream  of  the  base, 
with  full  compression  at  z/r  =  +2.0.   There  is  some  error 
in  the  calculation  introduced  by  neglecting  vorticity  and 
by  the  compression  fan  at  the  base,  rather  than  the  oblique 
shock.   See  Shaprio  [Ref.  31]  for  a  discussion  of  compres- 
sion fans.   Assuming  as  much  as  a  1°  error  in  the  flow  angle, 
far  exceeding  the  estimate  given  in  Section  IV,  the  value  of 
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Ph/Poo  remains  at  1.8  and  gives  specific  impulse  of  230 
seconds.   With  w/r  =  4.0  and  z/r  =  2.0,  part  of  the  com- 
pression due  to  combustion  was  applied  along  the  shear 
layer  downstream  of  the  critical  point.   With  z/r  =  1.0 
more  of  the  compression  was  applied  upstream  of  the  critical 
point,  but  the  base  pressure  was  less.   The  variation  of 
Vu/Voo   with  z/r  indicates  that  the  location  of  dp/ds  (due 
to  the  external  compression  waves)  along  the  shear  layer 
may  be  an  important  consideration  in  achieving  high  values 
of  base  pressure  with  external  burning.   This  hypothesis 
and  the  condition  where  the  external  flow  with  heat  addi- 
tion becomes  subsonic  require  additional  study. 

Figure  23  provides  some  details  regarding  the  optimum 
flow  condition  shown  in  Fig.  22  as  point  A.   The  heat  zone 
and  shear  layer  have  been  plotted  in  the  lower  half  of 
Fig.  23.   The  boundary  layer  thickness  has  been  included 
as  well  as  the  oblique  shock  wave  which  forms  when  p,  /p^ 
is  greater  than  one.   For  comparison,  the  shear  layer 
without  heat  addition  has  been  shown  as  a  dashed  curve  in 
Fig.  23.   Pressure  ratios  along  the  shear  layer  have  been 
plotted  in  the  upper  half  of  Fig.  23;  again,  the  dashed 
curve  is  without  heat  addition.   In  both  cases,  with  and 
without  heat,  the  shear  layer  pressure  ratio  shows  a  gradual 
increase  with  distance  downstream  from  the  base.   The  pres- 
sure ratio  along  the  inner  streamline  of  the  heat  zone 
(points  a1  through  e')  has  also  been  plotted  along  with  the 
shear  layer  pressure  ratios  (points  a  through  f ) .   This 
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streamline  experiences  an  abrupt  increase  in  pressure  ratio 
where  it  intersects  the  oblique  shock  wave  at  c';  then,  its 
pressure  ratio  starts  to  decrease  with  a  corresponding  in- 
crease in  Mach  number,  indicating  that  the  flow  along  the 
inner  streamline  accelerates  downstream  of  the  oblique 
shock  wave.   The  calculation  along  the  shear  layer  was 
terminated  before  reaching  a  maximum  pressure  ratio  due  to 
the  saddle  point  singularity,  Fig.  4,  since  the  straddle 
method  employed  in  the  computer  program  had  established  the 
required  base  pressure  ratio  within  ±  %% . 

Additional  calculations  with  the  computer  program  given 
in  Appendix  J  would  improve  the  accuracy  of  p,  /pOT  and  would 
take  values  for  the  shear  layer  slightly  closer  to  the 
critical  point,  but,  since  the  critical  point  is  a  saddle 
point  singularity,  would  not  be  expected  to  calculate  through 
the  critical  point.   Calculation  "through"  the  critical 
point,  using  the  straddle  method,  would  require  changes  in 
the  program  to  restrict  the  values  allowed  for  dF/dw   and 
force  the  calculation  across  the  singularity.   Additional 
study  is  required  to  include  this  change  and  provide  the 
capability  of  calculating  the  entire  base  flow  field. 

To  verify  the  existence  of  an  oblique  shock  wave  at 
the  base,  which  was  calculated  in  Fig.  23,  one  may  consider 
the  field  of  compression  waves  impinging  on   the  shear  layer 
as  being  concentrated  into  a  single  oblique  shock  wave  im- 
pinging on  the  shear  layer.   This  simplified  condition, 
shown  in  Fig.  24,  has  been  studied  experimentally  by  Korst 
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[Ref.*  36]  and  was  found  to  create  an  oblique  shock  wave  at 
the  base.   Thus,  the  calculated  shear  layer  profile  and 
pressure  ratios  in  Fig.  23  are  in  agreement  with  experimen 
tal  results  from  the  simplified  condition. 
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VI.   CONCLUSIONS 

With  the  method  of  characteristics  and  the  Crocco-Lees 
base  flow  model,  one  can  accurately  predict  the  base  pres- 
sure ratio  of  a  two  dimensional  or  axisymmetric  body  with 
moderate  external  burning  for  various  Mach  and  Reynolds 
numbers.   The  correct  values  for  the  two  transition 
Reynolds  numbers  are  essential  for  accurate  base  pressure 
ratio  calculation  of  a  particular  body  and  should  be  ob- 
tained from  experiment.   The  effect  of  projectile  spin  can 
be  included  approximately  by  dividing  the  transition 
Reynolds  numbers  with  one  plus  the  spin  number. 

Net  thrust  can  be  produced  by  external  burning.   Tlac- 
ing  the  combustion  zone  well  forward  of  the  base  to  give 
full  compression  on  the  body  was  found  by  experiment  to 
produce  a  decrease  in  thrust.   For  maximum  projectile 
thrust  with  supersonic  external  flow  the  computed  heat  zone 
placement  was  found  to  start  1.2  base  radii  upstream  of  the 
base  and  extend  4.0  base  radii  downstream.   The  zone  started 
with  an  inner  and  outer  streamline  of  2.0  and  3.0  base  radii 
respectively  and  received  a  fuel  flow  of  0.6  pounds  per 
second  of  fuel  rich  propellant  exhaust,  ARC- 168.   It  ap- 
peared higher  values  of  base  pressure  would  be  possible  with 
the  same  fuel  flow  rate  by  decreasing  the  length  of  the  com- 
bustion zone  still  further.   The  flow  field  could  not  be 
computed  with  this  heat  distribution  because  of  the  resulting 
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subsonic  flow  and  requires  additional  study  to  determine  any 
further  advantage  or  limits  under  these  conditions.   Experi- 
ments with  actual  combustion  should  be  conducted  to  deter- 
mine the  feasibility  of  the  desired  combustion  placement 
and  to  measure  the  resulting  base  pressure. 

The  rate  of  pressure  increase  along  the  shear  layer  due 
to  the  heat  zone  placement  appeared  to  influence  the  com- 
puted value  of  base  pressure  ratio.   This  effect  requires 
additional  study  to  determine  its  significance  in  achieving 
higher  values  of  base  pressure  ratio  with  external  burning. 
Experimental  studies  are  needed  to  determine  any  variation 
of  the  empirical  F  -  k  relation  when  external  pressure  waves 
impinge  on  the  shear  layer. 
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APPENDIX  A:   DERIVATION  OF  THE  CHARACTERISTIC  EQUATIONS  WITH 
HEAT  ADDITION 


Equations  (1)  through  (4)  govern  steady  compressible  flow 
with  heat  addition  [Ref .  18] .   These  equations  neglect  heat 
conduction,  viscosity  and  body  forces.   In  natural  coordi- 
nates, Pig.  1, U  =  ws .   The  continuity  equation,  Eqn.  (1), 
may  be  written 

(2iTr)   pwAn  =  constant  (1A) 

where  a  =  0  for  two  dimensional  planar  flow  and  a  =  1  for 
axisymmetric  flow.   The  momentum  equation,  Eqn.  (2),  may  be 
written  for  the  two  unit  vector  directions,  s  and  n, 

(  pw  |f  ♦  If  )s  =  .0  (2A) 

<  pw2  H  +  If  >»  ■  °  ■  <3A) 

The  energy  equation,  Eqn.  (3),  becomes 

3h 
~9s 


w*|  (h  +  i  w2)  =  q  =  w  S  (4A) 


where  the  total  enthalpy,  h~  =  h  +  yw2 ,  has  been  used.   The 

thermodynamic  relation  for  entropy,  S,  may  be  written 

TVS  =  Vh  -  -  Vp  (5A) 

9s      l    3s    p  dsJ            3s  K      J 

T9S      ,  3w    1  9p>.      0  nn 

ITT—  =   -    W7-  +   —  T^l   +   — Tv—   .  (7A) 

9n      v  3n    p  3i"r     9n  v  J 
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Using  Eqn.  (3A)  into  (7A) 

w  9n    8  s      w7"  8n    w2   8n  *"   •* 

Forming  the  logarithmic  derivative  of  Eqn.  (1A)  in  the  s  di- 
rection 

«  sine  +  I  IP  +  I  3*  +  3i  =  0  (9A) 

r         p  8?    w  8s    8n  v   ^ 

where  the  relations  8r/8s  =  sinO  and  [1/An]  [8An/8s]  =  80/8n 

have  been  used.   Forming  the  logarithmic  derivative  of  Eqn. 
(4)  in  the  s  direction 

I  l£  =  1  iP  +  I  II  (10A) 

p  8s    p  8s    T  8s  '  k    J 

Substituting  Eqns .  (2A)  and  (10A)  into  (9A) 

a   .  0  .  e    1    w.8w    18T    80    n  nin 

—  sine  +  (--T^F)T_"T^r-+T_=0-  (HA) 

r         v  w    Rf  J    8s    T  8s    8n  K         J 

From  Eqn.  (4A) 

q      8T     8w  Min 

21  =  c  -r—   +  Wtt—  (12AJ 

w     p8s     8s  v    J 

where  the  enthalpy  has  been  written  in  terms  of  a  constant 
c  ,  h  =  c  T.   Substituting  Eqn.  (12A)  into  (HA) 

,    1           w    w   n  9w    8  0      q     as  in  6  n       /11A-, 

f  —  -  ttt  +  — *v    )  -v-  +  v~  "  V  +  =  0»  (l^A) 

v  w    RT    c  I  J    8s    8n   wc  1      r  v 

P  P 

Combining  terms,  and  using  the  Mach  angle  relation  cot2y  = 
M2  -  1,  Eqn.  (13A)  becomes 

cot2y   8w    90  _  as i n0      q 


w    .8s    8n     r     wc  T 

P 
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(14A) 


Thfc  Prandtl -Meyer  Angle,  v,  is  defined  by  the  relation 


dv  =  £oty  dw 
w 


(15A) 


Using  Eqn .  (15A)  into  (14A)  and  multiplying  by  tan  y 

-4-  6A) 


3v  3G         .  I  as  me 

T?—   -    tany   -r—  =    tany  / - 

8s  3n  It  wc 


' 


TJ 


Similarly,  Eqn.  (ISA)  into  (8A)  becomes 
3v    36 


tany 


t  as  +  i   o 


3n    3s      w2  3n    w2   3n 


(17A) 


Adding  Eqns .  (16A)  and  (17A) 


(  7r|  +  tany  ~  )  (v  -  6) 
v  3s        3n  M 


tany 


asinO 


r     wc  T 
P 


* 


w2^  (  3n 


3h^) 
3n 


(18A) 


Subtracting  Eqn.  (17A)  from  (16A) 


3  3 

(  ^   tany  ^-  ) (v  +  6)  =  tany 


3s 


jasin8     q   V, 


w  { 


3S  _  8h0l 
w 2  f  x  3  n     3  n  I 


(19A) 


Now,  in  the  characteristic  coordinate  system 


secy 


n    ) 

3n 


(  37  +  tanvi  91T  )(  > 


(20A) 


and 


secy 


3(   )  - 

3C 


<  3?  ■  taniJ  3K  )(   } 


(21A) 
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With  this,  Eqns.  (18A)  and  (19A)  become 

(22A) 

3(v  -  0)         TasinG       q   /    cosy  J^dS         dh0   L 

L  an    =  siny  (~ wT/T  j  •  ^-  J1^  ■  -Tn") 

and 

(2  3A) 

»(v  +  °)  =  sinu  /asine   .  __g)  +  cosy  /  9S  _  ^0_f 
9£      Siny  f   r      wc  T  J  vi*      I   l  dn  9nJ' 

Eqns.  (22A)  and  (23A)  are  the  characteristic  equations,  with 
heat  addition,  governing  the  two  dimensional  planar,  u  =  0, 
or  axisymmetric ,  a  =  1,  flow.   These  equations  are  dimensional, 
having  the  dimension  of  reciprocal  length. 
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APPENDIX  B:   DERIVATION  OF  CONDITION  FOR  IRROTATIONALITY  WITH 
HEAT  ADDITION 


The  equations  required  for  this  derivation  are  momentum 
with  no  viscosity  or  body  force, 

Dt  =  "  -  Vp  (IB) 

where  U  is  velocity,  p  is  density,  and  p  is  static  pressure; 
energy  with  heat  addition,  but  without  viscosity,  body  force, 
or  conduction, 

Dffi^MD.  q  +  i  ||  (2R) 

where  h  is  enthalpy,  and  q  is  time  rate  of  heat  addition  per 
unit  mass;  the  second  law  of  thermodynamics 

TVS  =  Vh  -  -  Vp  (3B) 

where  S  is  specific  entropy;  and  the  vector  identity 

U-VU  =  VpiU-U)  -  TJX(VXU).  (4B) 

Combining  Eqns .  (IB),  (3B)  and  (4B) 

~  +  V(|  U-U)  -  UX(VXU)  =  TVS  -   Vh  .  (5B) 

Introducing  the  definition  for  total  enthalpy,  h„  =  h  +  ^U'U, 
into  Eqn.  (5B)  ,  one  obtains  Croccos's  Theorem 

TVS  +  UX(VXU)  =  VhQ  +  |y   .  (6B) 

Let    VXU   =    Q   and    take    the    curl    of  Eqn.     (6B) , 

JH  =    VTXVS    +    ft-VU    -    JTV'U  (7B) 
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since  V'ft  =  V* (VXU)  =  0  and  the  curl  of  a  gradient  is  zero. 
Eqn.  (7B)  describes  the  total  rate  of  change  of  vorticity, 

a. 

Now  assume  that  fi  =  0  prior  to  and  during  heat  addition. 
Eqn.  (7B)  becomes 

^  =  VTXVS  =  0.  (8B) 

Then,  Eqn.  (8B)  is  the  condition  which  must  be  satisfied  dur- 
ing heat  addition  if  the  flow  is  to  remain  irrotational ; 
that  is,  the  temperature  gradient  must  be  parallel  to  the 
entropy  gradient.   Referring  to  Eqn.  (3B)  and  assuming 

h  =  c  T,  where  c   is  constant,  the  condition  for  irrotation- 
p  p 

ality  may  be  stated  in  terms  of  pressure;  that  is,  the  tem- 
perature gradient  must  be  parallel  to  the  pressure  gradient. 

In  the  case  of  steady  flow,  the  irrotat ionality  condi- 
tion may  be  given  in  terms  of  heat  per  unit  mass  which  has 
been  added,  Q.   In  Eqn.  (2B)  let  q  =  U • VQ ,  and,  using  the 
total  enthalpy  hfl  =  h  +  ^U'U, 

U-V(h0  -  Q)  =  0.  (9B) 

Then  along  any  streamline 

hQ  -  Q  =  hQ  (10B) 

CO 

where  it  is  assumed  that  the  upstream  conditions,  before  'neat 
is  added,  i.e.  Q  =  0,  are  uniform  [Ref.  19].   With  these  as- 
sumptions one  may  take  the  gradient  of  Eqn.  (10B) 

V(hQ  -  Q)  -  0  CUB) 
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since  hn   is  constant.   Using  the  steady  form  of  Iiqn.  (6BJ  , 

00 

with  ?2  =  0,  and  Eqn.  (HB) 

VS  =  |  VhQ  =  ~  VQ.  (12B) 

Also,  from  Eqn.  (11B) ,  one  may  write 

VQ  =  V(h  +  yj-U)  =  c  VT  +  V(^U'U).  (13B) 

Combining  Eqns  .  (12B)  and  (13B)  in  (8B) 

—  /vQ  -  V(%U«tT)  J  X  ^   VQ  =  0.  (14B) 

Since  (VQ)X(VQ)  =  0,  Eqn.  (14B)  reduces  to 

(VQ)XV(%U'fT)  =  0.  (15B) 

Eqn.  (15B)  is  the  condition  which  must  be  satisfied  to  pre- 
serve irrotationality  in  steady  flow  with  heat  addition;  that 
is,  the  gradient  of  heat  added  must  be  parallel  to  the  gra- 
dient of  kinetic  energy. 

If  the  gradient  of  heat  added  is  a  small  perturbation,  as 
is  the  case  of  linearized  flow,  then  the  irrotationality  con- 
dition in  Eqn.  (15B)  is  automatically  satisfied.   To  show  this 
let  Q  =  h^Q',  where  |Q'|<<1,  and  U  -  U^ (I  +  g  ),  where  U » -U ' 

CO 

<<1,  then,  substituting  in  Eqn.  (15B) , 

VQXVC^U-U)  =  hooVQ,XV(Uoo^•lJ,  +  ^U-D"')  (16B) 

is  of  second  and  third  order  perturbation. 

This  condition  may  be  applied  to  Eqns .  (24A)  and  (25A), 
Appendix  A,  and  neglect  the  vorticity  term 


rT9s    d_\l 

\     9n     9n  J  ' 


- 
w 
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thus-  simplifying  the  method  of  characteristics  calculation. 
First  of  all,  if  the  heat  zone  is  started  in  a  region  of 
uniform  flow  along  a  line  normal  to  the  streamlines,  and  if 
the  heat  is  added  uniformly  in  a  zone  bounded  by  two  stream- 
lines, then  the  condition  of  Eqn .  (15B)  will  be  approximately 
satisfied  within  the  zone  since  the  gradient  of  heat  added 
will  be  approximately  parallel  to  the  gradient  of  kinetic 
energy.   Second,  if  the  first  condition  is  not  satisfied 
(i.e.,  along  the  streamline  boundaries  of  the  heat  zone  where 
there  is  a  jump  in  Q  normal  to  the  streamline),  the  vorticity 
may  still  be  neglected  if  the  heat  added  is  a  small  perturba- 
tion, that  is 

q  =  U-VQ  =  hJT-  VQ1  .  (17B) 

Then  the  heat  term  used  in  Eqns .  (24A)  and  (25A)  must  be 


m =  ^-  <<  1/r,. 


/r,.  (18B) 


wc  T      wc  T         3s        b 
P         P 

This  condition  is  satisfied  for  the  heat  zone  boundary,  where 
r,  q/wc  T  is  taken  as  0.171  to  get  -  4^°  deflection  of  the 
inner  heat  zone  streamline. 
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APPENDIX  C:   BRIEF  DEVELOPMENT  OF  CROCCO-LEES  BASE  FLOW  FORMU- 
LATION [Rcfs.  1,8,13] 


The  mass  flux,  m,  in  the  base  flow  is  given  by  the  inte 
gral 


0 


=  J^   (27rr)a  pudr  (1C) 


where  6  is  the  viscous  flow  thickness,  r  is  radius,  p  is  den- 
sity, u  is  axial  component  of  velocity,  a  equals  zero  for  two 
dimensional  planar  flow  and  a   equals  one  for  axisymmetric 
flow.   The  momentum  flux  is 


/ 


a  ~„2 


I  =  J        (27Tr)"  pu2dr  .  (2C) 

cr 

For  the  base  flow,  where  there  is  no  wall  drag,  the  momentum 
equation  becomes 

dl      dm    f    r-NOt  »  dp  r in\ 

cE  =  ue  dx  "  ^^       6   A  (3C) 

where  u   is  the  external  axial  component  of  velocity,  dp/dx 
is  the  static  pressure  gradient  and,  assuming  small  angles, 

^  =  (2Tr6)a  p  u  (  ££  -  6   ).  (4C) 

dx    v    '   Ke  e ^  dx    e  '  K      ' 

A  major  point  of  the  Crocco-Lees  theory  is  that 

^  -  6   =  k  (5C) 

dx     e  K      J 

where  k  is  a  non  dimensional  mixing  coefficient  which  must  be 
specified.   The  value  for  k  is  approximately  0.0033  for  laminar 
mixing  and  0.03  for  turbulent  mixing. 
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4  A  mean  velocity  is  defined  by 
Uj  =  I/in  (6C) 

and  a  velocity  ratio  by 

1'  e  v  J 

The  mean  density,  p,  ,  is  defined  by  the  equation 

m  :_-  ( 7i  6 )  '  5  p .  u  1  .  ( 8  C ) 

It  is  assumed  that  the  stagnation  temperature  is  constant 
throughout  the  flow  and  that  the  lateral  pressure  gradient, 
dp/dr,  is  zero.   Knowing  the  pressure,  p  =  p  ,  and  p1 ,  a 
mean  temperature,  T..  ,  can  be  defined 

T1    =    Pe/P1R.  (9C) 

Introducing  the  displacement  thickness,  6*,  and  momentum 
thickness,  6**,  of  ordinary  boundary  layer  theory 


Mx    s  n6  r.       .  (IOC) 

6*(l  +  a)    „        P      (2r)a(l    -    -EL.  )dr    =    6^1  +  a)  m 

ry  P     U  "     FT 

0  e   e  p   u   7Ta 

e    e 

and 


(1+a)    =    S6    (2r)a  -£2L-   (i    -    Ji  ) 

<T  pelIe 


(11C) 


(l+o)        ,*(l+o0  I 


6^    "J     -    6 


2    a 

p    u    I 
e    e 


Noting  that  k  =  I/mu   and  using  Eqns .  (IOC)  and  (1]C),  one 

e 

has 
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6O-  +  °0  _  6*(l  +  cx)    JC**(l  +  a) 

-  6* 


k  =  2 _6  ^J> .  (12C) 


On  application  of  the  Stewartson  transformation  [Rcfs.  8,13] 
one  finds  6,  6*,  and  6**  in  Eqn .  (12C)  are  replaced  by  the 
incompressible  quantities  6.,  6*  and  6s?*.   Thus,  within  the 
restrictions  of  the  Stewartson  transformation,  any  compres- 
sible laminar  base  flow  can  be  related  to  a  corresponding 
incompressible  laminar  flow.   It  is  assumed  this  also  applies 
in  the  turbulent  case.   The  ratio  of  mean  temperature,  T- , 
to  stagnation  temperature,  T  ,  may  be  written 

T 

r"  =  f  "    2    wi  (13C) 

where 


f  =  AH-*)\    6*(l  +  a)  (14C) 
i         i 

Wh  =  u./a  ,  and  a.  is  the  stagnation  speed  of  sound.  It  has 

been  found  more  convenient  to  work  with  the  quantity  F  >  0, 

F  +  1  =  f/K2 .  (15C) 

Experimental  data  [Ref.  8]  has  shown  a  relationship  between 
F  and  k  which  can  be  expressed  by  the  polynomial 

k    =    1.0    -    0.4F    +    0.16F2    -    0.01185F3.  (16C) 

For    F   >    1.5,    k    =    0.72    and    d</dF   =    0.0. 
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From  Eqn.  (8C)  one  may  write  an  expression  for  mass  flux 


m 


ma 


Y06)a    6    pewe 


k[F   +    1 


Y"l  2l 


(17C) 


Using  this  and  Eqns .  (3C)  ,  (4C)  and  (SC)  one  may  form  a 
single  ordinary  differential  equation 


dF 

dw 


rY  +  l        2 
v    2         e 


i)d 


2        e^ 


Lyw; 


kF 


1    - 


Y-l       2 

2      We 


(18C) 


1 


K 


W    (1-^iv2) 
e v         2      ey 


3tl  w2   e 


I.       _^1 

C  K      (1      -      K)     J 


(dK/dF)    -    1 


K     (1     "     KJ         k 

Eqn.  (18C)  has  a  saddle  point  singularity,  [Refs.  1,8],  also 
called  the  critical  point,  which,  in  integrated  form,  deter- 
mines the  one  unique  solution  describing  the  given  base  flow 
condition.   In  two  dimensional  planar  flow  without  external 

heat  addition,  where  9   is  known  from  w  by  the  Prandtl- 

e  e  ' 

Meyer  relation,  one  may  start  at  the  critical  point,  where 
both  the  numerator  and  denominator  of  Eqn.  (18C)  are  zero, 
and  integrated  toward  the  base  [Ref.  8].   In  flow  with  ex- 
ternal heat  addition,  and  in  axisymmetric  flow,  6   is  not 
known  as  a  function  of  w  ;  thus,  in  these  cases,  one  must 
start  integration  of  Eqn.  (18C)  at  the  base  and  concurrently 
solve  the  external  flow  field,  by  the  method  of  character- 
istics or  linearized  theory,  for  the  values  of  6^  along  the 
shear  layer.   This  process  is  one  of  trial  and  error  describe,'. 
in  section  IV  as  the  straddle  method. 
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APPENDIX  D:   CALCULATION  OF  BOUNDARY  LAYER  MASS  FLUX  AT  THE 
BASE 


The  boundary  layer  is  assumed  to  be  adiabatic,  have  con 
stant  pressure  across  its  thickness,  6,  and  have  a  velocity 

1  /n 

profile  (u/u  )  =  (y/<5)  '  ,  where  u   is  the  external  flow 
velocity  and  n  is  an  integer  [Ref.  24] ;  1  <  n  <  7,  for 
laminar  to  turbulent  flow.   The  boundary  layer  mass  flux  at 
the  base  is  given  by  the  integral 


Orb  +  6     a 
m,  =   I     (2irr)  pudr  (ID) 

rb 

where  a  is  zero  for  two  dimensional  planar  flow,  or  one  for 

axisymmetric  flow,  p  is  density  and  r,  is  the  base  radius  or 

half-height.   Using  the  equation  of  state,  p  =  pRT,  and  the 

above  assumptions,  Eqn.  (ID)  becomes 


b  r~Yu    l/n 

p     P         (2-rrr)au  (-r-5-)    dr 

m   =  _|    I     ^_2 (2D) 

bR  u 2    r  -  r,  ~  ,  v  J 

T    e  (  b-,  2/n 

lT    "  2c   l   6   J 
P 


xb 

J 


where  p   is  the  external  pressure.   Defining  w   as  u  /a  , 
where  a   is  the  stagnation  speed  of  sound,  Eqn.  (2D)  becomes 


r,  +6 

r-r. 


mb    ■   mbas    ■    V6pewc(2„6)a      [  (-^jl/n^^drj       (JD) 


2b  1         Y-l      2  fJ    rb.2/n 

1    "    T"  WeC~6~"J 
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Th'en 


m. 


1        1  /  ru 

'  x1/n  rx+4)a 


Mb  " 


Y<5pewe(27r6) 


a 


-t 


&+-%r  dx 


B2X2/n 


where  x  =  (r  -  r,  )/6  and  B2  =  ^5—  w2 
v     bJ  2    e 


n 


If  a  =  0  and  x  =  z  ,  Eqn .  (4D)  becomes 


D  =  *J  t-^ 


dz 


2D    J1  J   1  -  B2z2  * 
0 

This  may  be  integrated.   For  n  odd,  Eqn.  (5D)  gives 

(n-l)/2 


I 


UODD     I  2Bn  L  1-BZ 


y 

j^!    [n+l-2j]B 


2j 


(4D) 


(5D) 


} 


where  the  summation  is  omitted  if  n  =  1.   For  n  even,  Eqn 
(5D)  gives 


1    ,  rl+B  . 
JEVN      I  2Bn  L  l    b 


n/2  \(7D) 

^   * -T-l 

A\   [n  +  1  -  2j]BzJj 


n 


If  a  =  1  and  x  =  z  ,  Eqn.  (4D)  becomes 


AXI 


;2n,  ru  />    n 

z   dz  b  I*    z 

1-B2z^  6  n  J   1-B2 

0  0 


dz 

72 


(8D) 


The  second  integral  in  Eqn.  (81"))  is  the  same  as  Eqn.  (5D)  , 
which  is  evaluated  in  Eqns .  (6D)  and  (7D) .   Integrating  the 
first  integral,  Eqn.  (81))  becomes 
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XAXI 


=  n  [2B^+1  ^  ^  "  J   [2(n-j)\l]BzV 


j  =  l 

rh  (9D) 

+  — -  I 

The  results  of  Eqns .  (6D)  ,  (7D)  and  (9D)  are  used  with  Eqn . 

(4D)  to  evaluate  m,    or  m,    .   This  boundary  layer  mass 

D2D      DAXI 
flux  at  the  base,  before  the  corner  expansion,  m,  ,  is  set 

equal  to  the  mass  flux  after  the  corner  expansion,  m/ ,  ex- 
pressed in  terms  of  the  Crocco-Lees  formulation,  Appendix  C. 


yfi 'p'w' 

b    =  Y6PeWeI2D  =  mb    =  Vl   2  (10D) 

°2n     e  e  LU        D?.n   *-»  rF'  +  i  -1=-^'  2i 


2D       *  w  *"     "2D    k'  [F'+1--V^v'2] 

z   e 

Taking  k '  =  0.72  and  the  value  of  w1  assumed  from  the  straddle 
&  e 

method,  one  solves  Eqn.  (10D)  for  the  value  of  F'  which  starts 

the  integration  along  the  shear  layer.   Similarly,  in  the 

axisymmetric  case, 

7Ty6  '  2p  '  w' 

m,  =    27Ty62p   w    IAVT    =   m/  = — = (11D) 

bAXI  e    e    AXI  bAXI         k'[F'+1-    ^w«2] 

1  2      e    J 
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APPENDIX  E:   MASS  FLUX  CALCULATION  ALONG  THE  SHEAR  LAYER 

From  Appendix  C,  Eqn.  (17C)  ,  the  mass  flux,  m,  may  be 

written 

Y(TT6)a  6  p  w 
m  =  e  e  (IE) 

k[Y+1-^±-   w2] 
L      2    eJ 

where  a  =  0  for  the  two  dimensional  planar  flow  and  a  =  1  for 

axisymmetric  flow.   One  may  evaluate  Eqn.  (IE)  at  two  adjacent 

points  along  the  shear  layer,  m.  and  m.+1,  and  solve  for  the 

distance  between  these  two  points,  dx ,  by  using  the  additional 

relations  [Ref.  1,8]. 

-J}!  =  6  +  k  (2E) 

and 

&  =  (27i6)a  p  w   a2  (  %$-   -    6  ).  (3E) 

dx    *•    J      Ke  e   s  v  dx     J  K      J 

Assuming  small  angles,  Eqn.  (2E)  states  that  the  rate  of 

change  of  the  shear  layer  height  or  radius,  dS/dx,  is  equal 

to  the  flow  angle  0  plus  some  constant,  k.   Then  the  rate  of 

change  of  mass  flux,  dm/dx,  given  by  Eqn.  (3E) ,  is  the  mass 

entrained  by  the  shear  layer  from  the  external  flow. 

To  compute  m.  .,,  one  increments  F.  with  AF-,  say  AF.  = 
1  l+l '  i  l*    7    l 

0.1F.  with  the  condition  AF  >  0.08, 

F.+1    =    F.    +    AF. .  (4E) 

l+l  li  v      J 


Then 


w  =    w        +    AF./(dF/dw    )  (5E) 

l  +  l  i 
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where  dF/dw   is  computed  from  Eqn.  (18C).   The  value  for  < .  + 1 
is  computed  from  an  empirical  relation  using  F-+1>  given  by 
Eqn.  (16C)  ,  with  k  .  +  -.  =  0.72  if  F-+1  >  1.5.   The  isentropic 
relation  between  static  and  total  pressure,  in  terms  of  y 
and  w  ,  is  used  to  evaluate  p   at  i  and  i+1.   Average  values 
between  i  and  i+1  are  used  to  evaluate  d6/dx  and  dm/dx,  whore 
the  isentropic  relation  between  static  and  total  density  must 
be  used  with  the  latter.   Since 

m.+1  =  m.  +  (dm/dx)dx  (4E) 

and 

6i+l    =    6-    +    (d<S/dx)dx    .  (5E) 

One  may  solve  directly  for  the  value  of  dx,  in  the  two  dimen- 
sional planar  case,  which  corresponds  to  the  increment  AF- 

j  r         el       l  +  l-v/z-dO     dm  v  rar\ 

dxi,n  =  6i  <  sr  •  «— 7  )/(  as? "  as  >■  (6E) 

2D  l     l+l 

For  the  axisymmetric  case 

dx      =  \    (  dx     -  /  «   dx2     ).  (7E) 

XAXI    Z      1W        °i   dX     nAXI 

Eqn.  (7E)  may  be  iterated  several  times,  as  long  as  dx. 

XAXI 
<  1,  to  determine  the  correct  value  of  dx. 

XAXI 
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APPENDIX  F:   CALCULATION  OF  BOUNDARY.  LAYER  THICKNESS  AND 
MIXING  COEFFICIENT 


The  laminar  boundary  layer  thickness,  6.   ,  was  computed 
at  point  x  from  the  compressible  Blasius  solution  [Rcf.  8] 

6L  =  (5.0  +  1.2(Y-l)M2)x/(Rex)^2  (IF) 

where  x  is  the  distance  along  the  surface  from  the  leading 
edge  and  Re   is  the  Reynolds  number  based  on  this  distance. 
If  the  base  Reynolds  number,  Re,  ,  is  less  than  the  transi- 
tion Reynolds  number,  Re   ,  then  the  boundary  layer  at  the 
base  is  completely  laminar  and  its  thickness  was  computed 
from  Eqn.  (IF),  with  x  taken  as  the  body  length. 

If  the  base  Reynolds  number,  Re,  ,  is  greater  than  the 
transition  Reynolds  number,  Re   ,  the  boundary  layer  will 
be  partly  or  completely  turbulent.   The  turbulent  thickness, 
6T,  was  computed  from  [Ref.  8] 

6T  "  ■^T7T7!CcWcfiH«/«**)  w 

Cr.  is  von  Karman's  incompressible  turbulent  skin  friction 
coefficient  given  by  [Ref.  24] 

(0.242)/(C£i)^  =  log]0(RcxCf.)  C3F] 

and  C,.,  is  van  Driest 's  compressible  turbulent  skin  friction 
coefficient  given  by  [Ref.  24] 


?  U 

Q  1)  u  ~\     2 


(4F) 


iffi^V    sin"lB   =    logl0(RexC£M^    +    1.261og10(l-B2) 

(B    CfM) 
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7  2 

where  B   is  (y-l)w  /2.      Equations  (3F)  and  (4F)  must  be 

solved  by  iteration.   The  quantity  6/6**    (Appendix  C)  may 

be  evaluated  for  the  two  dimensional  1/7th  profile  using 

the  same  assumptions  employed  in  Appendix  D.   Briefly 

6**   "  a         T^T   C    1    -   ^-  )dr  [SF] 

r,  ^    e    eJ  e 


—j—  =    (1-— ~— w   )      \  k    6      J  v    6      J  f,    /Jf. 

6  2      G       O  "    y-1    2rr-r„2/7,         (dr/6) 

rb  (1      2   we("6-b)         } 


",2lf    1      mfl      1     -     V 


=    7  (1-*  )  (.-*■*  In(-^r) 


6  '  "2B8         Vb2  .^       (8-2j)B2J 


4 

~VnO    +     2  " 2T}-  (7F) 

2By         1    13  aj  (9-2j)BZJ 

j  =  1  v         J J 


Since  the  actual  transition  process  is  more  or  less  probabil 

istic  in  nature,  it  was  assumed  that  the  transition  from 

laminar  to  turbulent  flow  can  be  described  by  the  integrated 

normal  distribution.   It  was  necessary  to  assume  a  turbulent 

transition  Reynolds  number,  Re    ,  for  this  purpose,  where 

Re.  .  is  approximately  four  times  Re^  .   The  transition 
trt     lt  '  tr 

computation  was  made  with  the  error  function,  erf,  which  is 
exactly  related  to  the  integrated  normal  distribution.   The 
argument  of  the  error  function  was  formed  in  terms  of  base 
distance  from  the  initial  transition  point,  and  was  weighted 
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so  its  value  was  -0.96  at  Re.   and  0.90  at  Re.  . ,      A  numeri- 
4  tr  trt 

cal  span,  S,  of  4.0  takes  the  integrated  normal  distribution 
from  0.02  to  0.98.  Formulating  these  conditions  in  terms  of 
body  radius,  r, 

Ztr  =  (Rctr/Reb  "  1-0) (L/r)  (8F) 

is  the  transition  point  on  the  body  measured  from  the  base, 
and 

Ztrt  "  tRetrt/Reb  -  1-0)  &/r]  (9F) 

is  the  turbulent  transition  point.   Note  that  Z.    will  lie 

r  trt 

downstream  of  the  base  if  Re,.  .  >  Re, .   Forming  the  average, 

trt     b         ° 

Z«.     =  k{1<.      +  Z+  J  (10F) 

trav    2V  tr     trt-'  v  .  * 

is  the  distance  from  the  base  to  the  center  of  the  transi- 
tion zone.   If  Z.     is  zero,  the  integral  of  the  cumulative 

trav        '         & 

distribution  must  be  0.5.   The  value  of  erf  at  this  point 
will  be  0.0,  and  the  cumulative  distribution  integral,  Irn> 
is 

ICD  =  fc(l  +  erf).  (11F) 

The  values  of  erf  at  the  base,  ERF,  ,  range  from  -  1.0  to 

+  1.0.   In  terms  of  Z 

S     Z. 
ERF,  =  erf  (- -— (12F) 

b       2^  (Z.  ¥    -    Z„„  ) 

K    trt     tr  J 

Then,  when  Re,  >  Re^  ,  6T, . ,,-  is  given  by 
'         b      tr'   BASE     b  ; 

6BASE  "  WlU  +  ERFb)  +  WT(1  -  ERFb).  (15!') 
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In *Eqn.  (13F)  ,  ERF,  will  be  positive  if  Re,  is  less  than 

^(Re^.   +  Re^  4),  thus  as  Re,  ranges  from  Re.   to  Re.  .  ,  the 
ZK      tr     trt7  '  b     b  tr      trt ' 

contribution  of  5.  will  range  from  0.99  to  0.01  and  the  con- 
tribution of  6T  will  range  from  0.01  to  0.99.   The  value  for 
n  used  to  compute  the  velocity  profile  at  the  base  (see  Ap- 
pendix D)  is  formed  exactly  as  Eqn.  (13F) 

"BASE  "  «*  +  ERV  +  I  ^  "  ERV  (14F> 

since  n  ranges  from  1  to  7 . 

When  the  transition  zone  lies  along  the  shear  layer, 
one  must  compute  a  transition  value  for  the  mixing  coeffi- 
cient, k,  at  the  point  being  considered,  xRt  .   The  distance 
xRt  is  measured  from  the  base  in  terms  of  r.   To  carry  out 
this  calculation,  the  value  of  erf  at  xR,  ,  ERFR1  ,  is  written, 
as  in  Eqn.  (12F),  in  terms  of  Z 

S   (Ztrav  "  XBL} 
ERFm  =  erf(—  — ^^ —  .  (15F) 

v  trt     trJ 

Note  that  Eqn.  (15F)  reduces  to  (12F)  since  xRJ  is  zero  at 
the  base.   Using  Eqn.  (15F),  the  value  for  the  coefficient, 
k^y  ,  is  formed  in  a  similar  manner  from  the  laminar  and 
turbulent  mixing  coefficients,  ky  and  kT , 

kBL  =  likL(1  +  ERFBL)  +  ^M1  "  ERFBL)>  (16F) 

In  Eqn.  (16F)  ,  the  value  of  ERF~T  will  be  positive  if  Z 
is  downstream  of  x„, ,  and  vice  versa. 
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APPENDIX  G:   DERIVATION  OF  LINEARIZED  STEADY  PLOW  EQUATION 
WITH  HEAT  ADDITION 


The  equations,  (1)  -  (4),  are  required  for  this  deriva- 
tion. 

V- (pU)  =  0  (1G) 

pU-ViT  =  -  Vp  (2G) 

U-V(h  +  ^U-U)  =  q  (3G) 

p  =  pRT  (4G) 

Forming  the  directional  derivative,  U«V,  of  Eqn .  (4G) 

tJ-Vp  =  pRU'VT  +  RTU-Vp.  (5G) 

Substituting  Eqn,  (1.G)  into  (SG) 

pRU-VT  =  U-Vp  +  pRTV-U  (6G) 

Substituting  U'  of  Eqn.  (2G)  in  (6G) 

RU»VT  =  -   U-V(^U-U)  +  RTV-U  (7G) 

Assuming  c   is  constant,  using  h  =  c  T  and  c   =  yR/(Y-l), 
substitute  y/(Y"1)  times  Eqn.  (7G)  into  (3G) 


_1   - 

a' 


V-U  -  -i-  U«V(JsU-U)  =  q/c  T.  (8G) 

L  p 


Let  T  =  T  (1   +  T')  ,  q  =  h  q'   U  =  U  T  +  V<f>. 

With  this,  Eqn.  (8G)  becomes 

(UroT    +    V4>)-V 

V2cf>    -    -j (Wl    +    UJ>       +    JiV4>-V<|>)    -    q'/(l    +    T') 

af    (1    +    T') 
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Retaining  only  first  order  terms,  Eqn.  (9G)  becomes 


V2cj>    -    M2    <f>         =    q'  . 

r  oo     rxx  n 


(10G) 


This  is  the  linearized  steady  flow  equation,  with  heat  addi 
tion  for  irrotation  subsonic  or  supersonic  flow. 
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APPENDIX  II:   DETERMINATION  OF  THE  MACH  ANGLE 

This  appendix  derives  a  numerical  procedure  for  deter- 
mining the  Mach  angle,  y,  for  a  given  value  of  the  Prandtl- 
Meyer  angle,  v 

The  Prandtl -Meyer  angle  may  be  written  [Ref.  18] 


1  _   l,_,wa  *,k, 
g 


v  =  i  tan"x(g(M2-l)^)  -  tan_1(M2-l)  2  (1H) 


where  M  is  the  Mach  Number,  1  <  M  <  °°, 

g  -  C(Y-1)/(Y+D)%  (2H) 

and  y  is  the  specific  heat  ratio.   Let 

y  +  g  =  tt/2  (3H) 

where  y  is  the  Mach  angle  defined  by 

sin  y  =  1/M.  (411) 

From  Eqns.  (3H)  and  (4H) 

(M2-l)^  -  tane.  (5H) 

Using  equation  (5H) ,  equation  (1H)  may  be  written 

tang(v+c)  =  gtane.  (6H) 

In  Eqn.  (6H) ,  the  relation  between  v  and  e  is  still  implicit 
For  y  =  5/3,  g  =  h   and  Eqn.  (611)  can  be  reduced  to  the  ex- 
plicit form 

tan3  y  =  tany  .  (7H) 

Thus,  using  Eqns.  (311)  and  (711),  the  simple  relation  may  be 
written  when  y  =  5/3 
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H  =  |  -  2tan~1((tany)1/3).  (811) 

This  analytical  result  [Ref.  37]  was  first  given  in  terms 
of  M,  rather  than  u,  along  with  a  more  complicated  expres- 
sion when  y  =  5/4.   There  are  no  analytical  expressions, 
P  =  y(v,y)>  for  other  values  of  y>  consequently,  u  must  be 
determined  from  tables,  or  numerically,  in  all  other  cases. 

Equation  (6H)  may  be  written  in  a  form  for  iteration 
which  always  converges  to  the  correct  value  of  e 

ei+1  =  tan'1(|tang(v+ei))  (9H) 

for  0  <  e.  <  tt/2  and  0  <  v  <  v    .   The  value  v     is  given 
i  -  -    -   max  max     6 

by  Eqn.  (1H)  for  M  infinite 

Vmax  =  \   C(Cy+1)/(Y-1))%  "  1)  -  (10H) 

One  may  start  the  iteration  of  Eqn.  (9H)  with  an  approxima- 
tion for  e-  given  by  Eqn.  (7H) 

e1   =    2tan"1((tan^)1/3) .  (1111) 

This  approximation  fails  for  large  values  of  v    since  it  may 
give  a  value  e..  >  tt/2.   When  this  occurs,  set  £..  =  tt/2  to 
start  the  iteration  with  Eqn.  (9H) . 

While  Eqn.  (911)  will  always  converge  to  the  correct  value 
of  e ,  its  convergence  is  slow  and,  after  two  or  three  itera- 
tion, should  be  replaced  by  a  method  utilizing  the  derivatives 
In  effect,  Eqn.  (911)  is  used  to  find  an  approximation  of  e 
which  lies  in  the  range  of  convergence  of  a  faster  method. 
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One  method  which  provides  rapid  convergence  may  be  written 
from  Eqn.  (6H)  by  using  a  Taylor's  series  expansion  for  the 
separate  sides  about  z.    and  solving  for  Ac.   Let 

f1     =  tangCv+e^  (1211) 

i 

f2  =  gtane.  (13H) 

i 


then 


Ae2  Ae? 

Aci+i =  (fi.  -  £2.  +  -rf-  (£i:  •  **:>  +  -3} 


1     1 


(fill   -   f...)  +  -_-^£.     .   £•   ) 


(14H) 


where  Ae1  =  0  starts  the  iteration  in  Eqn .  (14H) .   The  advan- 
tage of  this  form  lies  in  the  fact  that  the  derivatives  may 

be  written  in  terms  of  the  values  of  f.   and  f~  ,  i.e.  f' 

i       i         i 
and  f'   involve  the  sec2x,  which  can  be  written  as  1  +  tan2x, 

i 
etc.   Thus  the  derivatives  may  be  included  with  only  two 

numerical  evaluations  by  the  computer  of  the  tangent  functions 
in  Eqns .  (12H)  and  (13H) .   Depending  on  the  number  of  deri- 
vatives included,  two  or  three  iterations  with  Eqn.  (14H) 
should  be  sufficient.   Then,  of  course, 

e..-  =  e.  +  Ac.  ,  (15H) 

l  +  l    1     j  +1  v    * 

is  used  to  evaluate  f1     and  f?    ,  etc.   For  y  =  1.4  and 

i+l       1  +  1 
including  the  fourth  derivative,  the  above  procedure  will 

provide  twelve  place  accuracy  for  e  in  two  iterations.   With 

the  required  value  for  c,  u  can  be  written  using  Eqn.  (311). 

This  method  was  used  in  the  computer  program  discussed  in 

Appendix  J. 
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APPENDIX  I:   SPINNING  CYLINDER  DESIGN 

The  spinning  cylinder  was  designed  to  be  compatible  with 
the  coaxial  free  jet  fixtures,  nozzles  and  spacers  used  for 
the  nonspinning  investigation  [Ref.  29].   The  spinning  cyl- 
inder and  its  support  assembly  had  the  same  external  dimen- 
sions as  the  stationary  cylinder  design,  1.282-inch  diameter 
and  9.95-inch  length,  to  insure  similar  test  results  without 
spin.   The  cylinder  was  driven  by  an  internal  impulse  tur- 
bine.  The  design  speed  of  48,000  revolutions  per  minute,  to 
match  the  spin  number  of  0.1595,  was  not  achieved.   This  was 
due  in  part  to  physical  interference  during  the  fabrication 
process  which  reduced  the  tangential  component  of  the  turbine 
jet  to  less  than  one-half  the  desired  value.   The  maximum 
speed  achieved,  without  nozzle  flow,  was  18,000  revolutions 
per  minute. 

The  spinning  cylinder  was  made  of  Plexiglas,  to  reduce 
rotating  mass,  and  was  glued  to  a  metal  endcap  which  formed 
the  cylinder  base.   The  endcap  was  driven  through  a  shaft 
which  screwed  into  the  internal  turbine.   The  shaft  also 
carried  thrust  loads,  from  an  internal  air  bearing  surface 
on  the  turbine  face,  to  react  with  pressure  and  friction 
forces  due  to  the  nozzle  flow.   This  part  of  the  design  was 
deficient  in  that  the  loads  due  to  nozzle  flow  caused  the 
turbine  face  to  rub  and  stopped  the  cylinder.   Time  did  not 
permit  modification  to  correct  this  deficiency. 
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The  spinning  cylinder  was  positioned  radially  by  24  air 
bearing  pads  provided  along  the  outer  surface  of  the  support 
cylinder.   High  pressure  air  was  supplied  to  these  pads 
through  orifices  formed  by  drilling  into  supply  tubes  sol- 
dered into  the  walls  of  the  support  cylinder,  Fig.  11.   The 
bearing  air  exhausted  through  ports  to  the  inside  of  the 
support  cylinder  where  it  mixed  with  turbine  exhaust  and 
flowed  out  the  bottom. 

The  base  pressure  on  the  spinning  cylinder  was  measured 
by  a  total  pressure  probe,  positioned  on  centerline  next  to 
the  base,  and  a  mercury  manometer.   The  cylinder  speed  was 
determined  through  a  photo  diode  and  a  pulse  counter.   The 
photo  diode  was  illuminated  by  a  laser  beam  reflecting  from 
the  metal  endcap  and  half  of  the  endcap  was  painted  black  to 
generate  one  pulse  per  revolution. 


Figure  II.   Disassembled  View  of  Spinning  Cylinder 
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APPENDIX  J:   COMPUTER  PROGRAM 

The  FORTRAN  computer  program,  included  at  the  end  of 
this  appendix,  required  150K  for  compile,  110K  for  link  and 
80K  for  execution.   It  was  written  to  compute  base  pressure 
ratio,  within  approximately  O.S  percent,  by  applying  the 
equations  and  methods  discussed  in  this  thesis.   In  view  of 
the  empirical  information  and  approximations  employed  in 
many  of  the  equations,  computation  to  greater  accuracy  did 
not  appear  warranted.   Using  an  IBM  360  computer,  computation 
of  the  base  pressure  ratio  for  one  flight  condition  required 
approximately  two  minutes  for  the  fully  turbulent  shear 
layer  without  heat  addition.   The  addition  of  heat  upstream 
of  the  base  and  low  Reynolds  number,  which  involved  a  long 
laminar  shear  layer,  increased  the  size  of  the  characteristic 
mesh  and  resulted  in  computation  times  as  high  as  ten  minutes 

The  flight  condition  was  specified  before  making  a  com- 
puter run.   This  included  choosing  the  Mach  number,  card  80, 
and  the  base  Reynolds  number,  card  140.   The  two  no  spin 
transition  Reynolds  numbers,  cards  120  and  130,  were  choosen 
for  the  particular  body  to  provide  agreement  with  available 
experimental  data.   The  ratio  of  body  radius  to  body  length 
was  specified  on  card  150.   The  spin  number,  rw/U,  was  given 
on  card  131.   The  ratio  of  specific  heats  for  air  was  set  at 
1.4  on  card  200.   Card  ISO  and  190  were  interchanged  for  two 
dimensional  planar  flow.   For  flow  without  heat  addition, 
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card  90  was  set  to  zero,  card  60  to  1STRT  =  90,  and  card  100 
to  ZSTRT  =  20.0. 

Several  additional  considerations  were  required  to  com- 
pute the  flight  condition  with  external  burning.   The  start- 
ing values  for  the  radii  of  the  heat  zone  inner  and  outer 
streamlines  were  specified  on  cards  420  and  430.   The  axial 
positions  to  start  and  stop  heat  addition  were  given  on 
cards  100  and  110,  with  reference  to  the  base.   The  value 
for  Q,  card  90,  was  specified  to  obtain  a  desired  -  4%°  de- 
flection of  the  inner  heat  zone  streamline,  however,  this 
value  could  have  been  specified  by  assuming  various  fuel 
flow"  rates  and  heating  values.   Finally,  the  value  for  ISTRT, 
card  60,  was  set  low  enough  (minimum  value  of  1)  to  place 
the  Mach  line  boundary  condition  well  upstream  of  the  outer 
compression  waves  generated  by  the  heat  zone.   The  values  of 
selected  quantities  were  printed,  to  monitor  the  progress  of 
the  computation,  in  addition  to  the  required  values  for  the 
base  pressure  ratio. 

A  simplified  flow  diagram  is  shown  in  Fig.  1 J .   It  is 
intended  to  show  the  basic  sections  of  the  computer  program 
and  computational  details  have  been  omitted. 

In  the  computer  program  the  group  of  IF  STATEMENTS  which 
examine  the  mesh  point  for  heat  addition  can  introduce  a 
slight  error  at  two  corners  of  the  heat  zone.   The  amount  of 
heat  assigned  between  two  points  is  proportional  to  the  ele- 
ment of  length  along  the  characteristic  in  the  heat  zone. 
The  IF  condition  is  checked  in  radial  position  against  the 
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START 


1 

ESTABLISH  PLIGHT  CONDITION  VARIABLES: 

M,  Reb,  Retr,  Retrt,  Y,  kL*  kT ,  r/L,  a,  Heat 

Zone  (Q,Z  .  . ,Z    +       ,R-  ,R   ),  Spin,  I  .   .  for 
^x'  start*  stop'  is'  os;'  *             start 

Characteristic  Mesh  Upstream  of  Base  to 

Enclose  Heat  Zone  (Z    ,   .   .  =  -0.1(100-1)). 

v  mesh  start        ^     *  J 

Compute  Base  Boundary  Layer:  r/6    and  n. 

Compute  Characteristic  Mesh  Along  Body  and 
Through  Heat  Zone  Up  to  the  Base.   Store  Values 
Along  Base  Characteristic  to  Reset  for  the 
Straddle  Method. 

Compute  Mass  Flux  at  the  Base. 

Assign  Value  of  w   for  Straddle  Method 

e 

Compute  Expansion  Fan  and  F. 


Compute  Increment  Along  Shear  Layer  from 
dF  =  -0.1F,  dw   =  dF/ (dF/dw  )  and 

m.  .  =  m.  +  (dm/dx)dx. 
l+l     l    y      i        j 


Compute  Characteristic  from  Shear  Layer  to 
Upstream  Mach  Line.   Correct  9  at  Shear  Layer 


Proceed  Along  Shear  Layer  Until  Critical  Point 
Where  dF/dw   Becomes  Negative.   Compute  p,  /p^ 

from  Value  of  w   at  Base.   Assign  Next  Value  of 

e  6 

w      at    Base    from   Sign    of   dF/dw      Numerator. 
e  &  e 


■STRADDLE    LOOP  <■ 


STOP 


Figure  1J.   Simplified  Plow  Diagram  for  Externa]  Burning  Assist 
Projectile  (EBAP)  Computer  Program, 
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Inner  and  outer  streamlines  of  the  heat  zone,  and  in  axial 
position  against  the  start  and  end  of  the  heat  zone.   Thus 
a  right  running  element  at  the  lower  left  or  upper  right 
corner  of  the  heat  zone  could  be  outside  the  zone,  but 
still  pass  the  IF  condition  tests  and  be  assigned  \   value 
for  heating.   This  error  is  small  with  a  small  characteristic 
mesh  and  is  neglected.   The  same  applies  to  a  left  running 
element  at  the  upper  left  or  lower  right  corner  of  the  heat 
zone . 
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